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For a five-parameter discrete 4 model, we derive various exact static solutions, including the 
staggered ones, in the form of the basic Jacobi elliptic functions sn, cn, and dn, and also in the form 
of their hyperbolic function limits such as kink (tanh) and single-humped pulse (sech) solutions. Such 
solutions are admitted by the considered model in seven cases, two of which have been discussed in 
the literature, and the remaining five cases are addressed here. We also obtain sine, staggered sine 
as well as a large number of short-periodic static solutions of the generalized 5-parameter model. All 
the Jacobi elliptic, hyperbolic and trigonometric function solutions (including the staggered ones) 
' are translationally invariant (TI), i.e., they can be shifted along the lattice by an arbitrary xo, but 

among the short-periodic solutions there are both TI and non-TI solutions. The stability of these 
• solutions is also investigated. Finally, the constructed Jacobi elliptic function solutions reveal four 

' new types of cubic nonlinearity with the TI property. 

00 '. PACS numbers: 05.45.-a, 05.45.Yv, 63.20.-e 

XJ} . I. INTRODUCTION AND SETUP 

d : 

In recent years there has been a growing interest in the analysis of new discrete nonlinear models since they play 
c~| \ a very important role in many physical applications. For example, the question of mobility of solitonic excitations in 
l— ~~ '■ discrete media is a key issue in many physical contexts; for example the mobility of dislocations, a kind of topological 
solitons, is of importance in the physics of plastic deformation of metals and other crystalline bodies Similar 
questions arise in optics for light pulses moving in optical waveguides or in photorefractive crystal lattices (see e.g., Q 
for a relevant recent discussion) and in atomic physics for Bose-Einstein condensates moving through optical lattice 
potentials (see e.g., Q for a recent review). These issues may prove critical in aspects related to the guidance and 
■ manipulation of coherent, nonlinear wavepackets in solid-state, atomic and optical physics applications. 

In particular, the translationally invariant (TI) discrete models [1] have received considerable attention since they 
admit static solutions that can be placed anywhere with respect to the lattice. Such discretizations have been 
— i | constructed and investigated for the Klein-Gordon field 0, H, IS 0> B B Gil EL [SElS EH, EH an d f° r the nonlinear 
Schrodinger equation 0, E3, EM EU H3, HH • For the Hamiltonian TI lattices [a^lll, this can be interpreted as the 
absence of the Peierls-Nabarro (PN) potential [![. For the non-Hamiltonian lattices, the height of the Peierls-Nabarro 
^ , barrier is path-dependent but there exists a continuous path along which the work required for a quasi-static shift of 
the solution along the lattice is zero [e|. 

In general, one can state that coherent structures in the TI models are not trapped by the lattice and they can 
be accelerated by even a weak external field. This particular property makes the TI discrete models potentially 
interesting for physical applications and one such physically meaningful model has been recently reported (l3| . 

For some of the TI models it has been demonstrated that they conserve momentum Q or energy (Hamiltonian) 
@, 3 (see also [H El)- 

However, we do not know a TI model conserving both momentum and Hamiltonian and, 
for the Klein-Gordon lattices with classically defined momentum, it was proved that these two conservation laws are 
mutually exclusive [10 ]. 

It may be noted that the TI discrete models can support even moving solutions, but only for selected propagation 
velocities 19]. In some cases, the exact static or even moving solutions to the TI models can be expressed explicitly in 
terms of the Jacobi elliptic functions (JEF). Even in the cases when JEF solutions are impossible, static solutions to 
a TI model can always be obtained iteratively from a nonlinear map (first integral) , solving at each step an algebraic 
equation. 

While there has been no universally acceptable definition of TI models, it is fairly easy to describe what a TI solution 
is. It is a static solution which can be placed anywhere with respect to the lattice. In particular, if there is an analytic 
TI solution with an arbitrary shift xq along the chain or if one can show that there is a corresponding Goldstone mode 
with zero frequency for any xq. As far as TI models are concerned, it is believed that they should possess following 
properties (i) they admit static solutions which can be placed anywhere with respect to the lattice, which can be 
associated with the absence of PN barrier. Note that in case analytical static solutions can be constructed with an 
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arbitrary shift along the chain, xq, that would automatically imply the absence of the PN barrier (ii) static version 
of TI discrete models are integrable, i.e. the static problems are reducible to a first-order difference equation which 
can be viewed as a nonlinear map from which static solutions can be constructed iteratively (in this study we will 
show that non-integrable three-point static problems also can have particular TI solutions derivable from a set of two 
lower-order finite-difference equations, and one of those equations is a two-point one while another is a three-point 
one) (iii) static solutions in TI models possess the translational Goldstone mode with zero frequency for any xq. 

A prototype class of discrete models, relevant to a variety of applications are the so called discrete <fi A models which 
feature a cubic nonlinearity. The purpose of this paper is to study in detail several issues related to TI models. 
In particular we consider a rather general discrete <fi A model with cubic nonlinearity which is invariant under the 
interchange of <f> n +i and </>„_i 

4>n = Jo (4>n+l + 4>n-l - 20n) + A0„ - A\<j>\ TT^ntyn+l + 4>n-l) 

-^-4>n{4>n+l + 4>l-l) - - ~^ ± ' ' 

with the model parameters satisfying the constraint 



n\Y n +l + 4>n-l) ~ A±<fi n <fi n+1 <fi n ^ X - — 0„+10„-l (</>n+l + 4>n-l) ~ "^"(^n+l + ( l ) 



E^ = A - ( 2 ) 



fc=l 

In Eq. JT]), <fi n (t) is the unknown function defined on the lattice x n — hn with the lattice spacing h > and overdot 
means derivative with respect to time t. Without the loss of generality it is sufficient to consider A = 1 or A = — 1. 

If model parameters Ak are constant (i.e. independent of h), condition of Eq. © ensures that, in the continuum 
limit, Eq. ([1} reduces to the <fi A equation 

4>tt = <t>xx + A-Xi - 4> 2 )- (3) 

On the other hand, if model parameters Ak are functions of h, then the continuum limit can be different from Eq. 
© even when Eq. © holds. 

Static form of Eq. ^ has the first integral 

<^-^(W 2 ) 2 + C = 0, (4) 

with the integration constant C. 

So far as we are aware of, all the discrete <fr models discussed in the literature, under static consideration, are 
special cases of the general model, Eq. (Q}. Some of these models are 

Model 1. Only A\ nonzero with other Ak being equal to zero results in the classical discretization of Eq. that 
has received a great deal of attention from the researchers in various fields. This model is not a TI one and it will not 
be further discussed here. 

Model 2. Ax = A 3 /2 — A± = A<5, A 5 = A 6 = 2Xj, A 2 = A(l — 47— 4<S) with arbitrary 7 and S. This non-Hamiltonian 
TI model (for arbitrary 7 and S) conserves momentum [4] 

P. 1 = 0n(0n+l - 0n-l ) • ( 5 ) 

Static version of this model (with the omitted inertia term (/>„) has the first integral with the integration constant C 

U((f> n -x, (j) n ) = ((/) n - <fi„-i) 2 + A(£ n n _i - A7 (0* + (fin-l) 

-M^n-i (<ft + 4>l-i) - A Q - 2 7 - 2/) <fil<fil_ x - ^ = , (6) 



with A defined by 



A = \h 2 , (7) 



from which any static solution to Eq. |T]) can be constructed iteratively, starting from any admissible value of (fio and 
solving at each step the algebraic problem. This is so because Eq. |T]) is nothing but 

- U(<p n ,<t> n+ i) - U(<f> n -i,(f) n ) 

"P™ = 1 1 ' ( 8 ) 

<Pn+l — <Pn-l 
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Equation ([6]) is the discretized first integral (DFI) [l2j|, i.e., in the continuum limit (h — > 0) it reduces to Eq. ((4]). DFI 
of Eq. ([5]) is quartic in both 4> n -\ and <^>„ and thus, it cannot be reduced to the integrable nonlinear map reported 
in (22J where the corresponding first integral is quadratic in </>„_i and 4> n - In this complete form the model was first 
constructed in Q and almost concurrently in 0] , although in the latter work any relation to the DFI was not observed. 
It is worth noting here that since 7 and 8 are arbitrary, one really has three distinct models each of which (as well as 
their sum) conserves momentum P\ and there is a corresponding two-point map in each case. 

Model 3. Model 2 with 7 = i5 = is the Bender-Tovbis model [|| . In the framework of the DFI approach [e| , 
almost the entire space of static solutions supported by this model was described and many of those solutions were 
expressed in terms of the Jacobi elliptic functions ( JEF) [l2[ ■ 

Model 4. Model 2 with S = and 7 = 1/4 is the Kevrekidis model [J. 

Model 5. Discrete </> 4 model 

$ n = ^2 (^"-i ~ 2( ^" + ^ > "+ 1 ) + V- A^ 2< /2 ' ^ 

discovered in does not belong to Eq. ([I]) but its static problem coincides with that of the Bender-Tovbis Model 3 
[1, HH • Some very special features of this TI model are the conservation of energy and the on-site discretization of 
the nonlinear term. In all other TI </> 4 models derived so far the nonlinear term is discretized on the three neighboring 
nodes (i.e. lattice sites). 

Model 6. With only A4 nonzero and other Ak being equal to zero, one arrives at the non-Hamiltonian TI model 
derived by Barashenkov, Oxtoby, and Pelinovsky Q and referred to as BOP. This model conserves the momentum 
defined as [lj] 

P2 = Y,4>n{<t>n+2-<l>n-2). (10) 
n 

The first integral of the static version of this model has been found in [l4| , where an almost complete set of static 
solutions supported by this model were derived and many of those solutions were expressed in terms of JEF. 

Model 7. Taking A\ = 2A/9, A 2 = A 3 = A/3, A± = A 5 = 0, A e = A/9 one gets the Hamiltonian of the Speight and 
Ward (SW) model @. For this model it is known that, it supports TI static kinks derivable from the two-point map 




6- A ' v ' 



where one can take either the upper or the lower sign. Note that this map is defined in case < A < 6. To get this 
map, one has to set in Eq. ^ C = and present it as 



6 a ±J^(l-0 2 ) = O. (12) 



Discretizing Eq. (fT2"|) as 



1 



H y/2 3\/2 



= 0, (13) 



and solving the resulting quadratic equation we come to Eq. (|11[) . It is not known if this model supports static TI 
solutions other than the kink. It is also not known if this model has the first integral of the static problem apart from 
the case of C — 0. In the present study we will give evidence that the answer to the second question is negative (see 
Sec. IVII Ap but we were able to find other TI solutions to this model, see Eq. (j9"7) and Eq. (j9"5j) . 

Model 8. A x = 4Aa( 7 +0),A 2 = 4A[2a 2 + 7 2 + /3( 7 - a)], A 3 = 4Aa(2 7 + /3), A 4 = 4A 7 (a -0),A B = 4Aa(a - fl), 
Aq = 4Aa 2 , with two free parameters a and (3 with 7 = 1/2 — 2a. This model was also offered by Barashenkov, 
Oxtoby, and Pelinovsky . This model includes as special cases the Model 3 (at a = (3 = 0), the Model 6 (at a = 
and (3 = -1/2) and also the Model 7 (at a = (3 = 1/6). 

Model 9. A x = 0, A 2 = 2A(l/2 - /3), A 3 = Acr(4 + h 2 ), A A = 2A/3 + Act(4 + h 2 ), A 5 = -8Acr, A 6 = 0, with two 
free parameters (3 and a. This model was also proposed by Barashenkov, Oxtoby, and Pelinovsky Q. Note that the 
coefficients A 3 and A4 in this model are /i-dependent and that the constraint Eq. ([2]) for this model is satisfied only 
in the continuum limit (i.e. for h = 0). 

In this paper we shall discuss four more TI models with cubic nonlinearities. The paper is organized as follows. 
In Sec. |TI]we discuss the subclasses of model Eq. ([IJ that support different conservation laws. In Sec. Mil we report 
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on a number of TI static solutions to Eq. ([I]) expressed in closed analytical form. All seven cases, when model Eq. 
(H|) supports the exact static JEF solutions, are described. Two of these seven cases have been previously studied 
in the literature, and for the remaining five cases, basic JEF solutions are given here together with their hyperbolic 
function limits. We also obtain a periodic sine solution for the general 5-parameter model as given by Eq. |T]). Section 
IIVI presents a number of short-periodic exact static solutions. In Sec. [V] we discuss the two-point maps for some of 
the TI models that have been reported in the literature and derive the map for the linear combination of Models 
3 (Bender- Tovbis) and 6 (BOP). Goldstone translational modes of the TI static solutions are discussed in Sec. I VII 
Numerical results that illustrate some important properties of the TI static solutions are presented in Sec. IVII1 The 
discussion, conclusions, and future challenges are described in Sec. IVIIII Finally in Appendix we spell out the short 
periodic solutions that are admitted by many of the models discussed in this paper. 



II. MOMENTUM AND ENERGY CONSERVATION 



As noted above, for the discrete model of Eq. Jl}, the momentum operator as given by Eq. © is conserved provided 
only A 2 is nonzero (Model 3), or if only A5 and A§ are nonzero with A5 — A§ (Model 4), or if only Ai,A3,A4 are 
nonzero with A\ — A3/2 ~ A4 and for a linear combination of these models which is precisely the Model 2. It is also 
conserved in Model 8 (provided a = f3 = 0) and in Model 9 (provided a = (3 = 0). On the other hand, as was already 
mentioned, the momentum defined by Eq. (| 1 Ojl is conserved for Model 6 with only A4 nonzero. It is also conserved 
in Model 8 (provided a = 0,/3 = -1/2) and in Model 9 (provided a = 0, /3 = 1/2). 

Model 10. For 



A x = 4aiA, A 2 — 6a 2 A, A 3 = 4a 3 A, A 6 = 2a 2 X, A 4 = A 5 = 0, with a\ + 2a 2 + a 3 = - , (14) 
the model Eq. ([1]) has the Hamiltonian 



{4>n 



,-l) 2 



2h 2 



A 
4 



A 
2 



an 



a 2 (, 



(15) 



and hence energy is conserved in this model. As it can be seen, the Hamiltonian model has two free parameters. Note 
that the Speight- Ward model 7 is a special case of this model for a\ — a 2 — 1/18 and a 2 — 1/12. Similarly, energy 
is also conserved in Model 8 provided a = /3 = 1/6. 



III. TRANSLATIONALLY INVARIANT JEF, HYPERBOLIC, AND TRIGONOMETRIC STATIC 

SOLUTIONS 

We shall first discuss the JEF solutions as well as the hyperbolic solutions which follow from the JEF solutions and 
later on we shall discuss the sine solutions which exist in almost all the models. 



A. JEF and Hyperbolic Solutions 

We shall now show that JEF solutions can be obtained for the discrete model of Eq. (JlJ in case A\ = Aq = in 
the following seven cases: (i) only A 2 nonzero; (ii) only A4 nonzero; (iii) only A 2 and A4 nonzero; (iv) only A3 and 
A§ nonzero; (v) A 2 , A3, and A5 nonzero; (vi) A3, A4, and A5 nonzero; (vii) A 2 , A3, A4, and A5 nonzero. 

The JEF solutions have already been reported in the case (i) in and in the case (ii) in [lj]. In this paper we 
report on the JEF solutions for the cases (iii) to (vii). In particular, in all these cases we obtain the sn solution 

4> n = SAsri[h/3(n + xq), m] , (16) 

the cn solution 

4> n = SAcn[h/3(n + xo), m] , (17) 

and the dn solution 



4> n = SAdn[h/3(n + xq), m] , 



(18) 
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where xo is an arbitrary shift, < m < 1 is the JEF modulus, and S — 1 for the non-staggered and S = (— l) n for 
the staggered solutions. 

We also derive the hyperbolic function limits (m = 1) of the above JEF solutions. In particular, while the sn 
solution reduces to the kink solution of the form 

(f)„ = SA tanh[h(3(n + x )], (19) 

both cn and dn solutions reduce to the single-humped pulse solution of the form 

(f>„ = SAsech[h(3(n + x )\ ■ (20) 

Expressions that relate the solution parameters A and [3 to the model parameters Aj- and, where applicable, the 
relations between Ak, are given in what follows. 

Note that if <j> n (t) is a solution to Eq. (TTJ), then the staggered solution (— l) n <f> n (t) is also a solution to the same 
equation, i.e. Eq. ([I]) with the coefficients A% , A3, and A4 having the opposite signs and further 2 — A is to be replaced 
by A — 2, where A is as given by Eq. ([7]). To make the presentation of the results as compact as possible, in most 
cases, we shall therefore not give the parameters for the staggered solutions. Only for the case (iii) with A2 and A4 
nonzero and for hyperbolic solutions, parameters will be given for both the staggered and the nonstaggered solutions. 

B. Case (iii): only A2 and A4 nonzero 

Since only A 2 ,A 4 are nonzero, hence Eq. reduces to 

A 2 + A 4 = X. (21) 

In this case, one has both nonstaggered sn and staggered sn solutions Eq. (fT6| . with S — 1 and S = (—1)™ respectively, 
provided the following relations are satisfied: 

2m 

= A 2 ns 2 (h/3,m) ± 2A 4 ns(hf3 , m)ns(2hp , m) , (22) 



A 2 h 2 

± (2 ^ 2 ^ 2 )m - A 2 cs(hf3)ds(hf3) ± A 4 ns 2 {hp, m) , (23) 

where ns(x,m) = l/sn(x,m), cs(x,m) = cn(x , m) / sn(x , m) and ds(x,m) = dn(x,m)/sn(x,m). In the above equa- 
tions, the upper sign corresponds to the nonstaggered sn solution (S = 1), while the lower sign corresponds to the 
staggered case (S — (— 1)™). 

In the limit of m = 1 the nonstaggered solution reduces to the kink solution, Eq. (|19p , with S = 1 and the relations 
(|22| and (|23|) take the simpler form 

A» = l, ^ 4 = 2tanh 2 (^)-A = _ 

tanh 2 (h/3) V ' 

Note that this solution is valid for any value of A including A = 2. 

On the other hand, in the limit of m — 1, the staggered sn solution reduces to the staggered kink solution, Eq. 
(fT9|) . with S = (-1)" and the relations j22]) and (J23J) reduce to 

Another set of exact JEF solutions are the nonstaggered as well as the staggered cn solutions as given by Eq. (fT7|) 
with S = 1 and S = (—1)", respectively, provided the following relations are satisfied: 

2m o 

— — = -A 2 ds 2 (h0, m) =f 2A 4 ds(/i/3, m)ds(2/i/3, m) , (26) 
A z h l 

± ^ 7 2 ^ m = T ^ds 2 [hP, m) - A 2 cs(ft,/3, m)ns(h(3, m) . (27) 
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In the above equations, the upper sign corresponds to the nonstaggered sn solution (5 = 1), while the lower sign 
corresponds to the staggered case (5 = (—1)™). 

Yet another exact JEF solutions are the nonstaggered as well as the staggered dn solutions as given by Eq. (fl8|) 
with 5 = 1 and 5 = (—1)" respectively provided the following relations are satisfied: 

o 

-A 2 cs 2 (/i/3, m) =f 2A*cs(fy3, m)cs(2h(3, m) , (28) 



A 2 h 2 

± ~ ^ = TAiCS 2 {hp, m) - A 2 ns(h/3, m)ds(h(3, m) . (29) 

In the above equations, the upper sign corresponds to the nonstaggered sn solution (5 — 1), while the lower sign 
corresponds to the staggered case (5 = (—1)"). 

In the limit of m — 1, the nonstaggered cn as well dn solutions reduce to the pulse solution Eq. ([20]) with 5 = 1 
and the relations (|26|) and ([27]) as well as ([28]) and ([29]) take the simpler form 

A-^W)-1]<0, 

Thus the pulse solution exists only if A < 0. 

On the other hand, in the limit of m = 1 the staggered cn as well dn solutions reduce to the staggered pulse solution 
Eq. ([20]) with 5 = (-1)™ and the relations {26]) and {27]) as well as {28]) and {29]) take the simpler form 

A = 2[cosh(^) + l]>0, A> = ^-^l y (31) 

Thus the staggered pulse solution ([3"T]) exists only if A > 4. 

C. Case (iv): only A3 and A$ nonzero 

Since only A 3 , A 5 are nonzero, Eq. {2]) reduces to 

A 3 + A 5 = X . (32) 
In this case, one has the sn solution Eq. {TB]l with 5 = 1 provided the following relations are satisfied: 

A 5 ns(2h(3, m) = -A 3 as(hl3, m) , (33) 

2m 

= -A 3 cs{hf3,m)ds(hf3,m) - A 5 [cs(2h(3,m)ds(2hp,m) - ns 2 {2hf3,m)} , (34) 



A 2 h 2 



{2 J^ m =-A 3 ns 2 (h(3,m). (35) 
In the limit m — 1, this reduces to the kink solution Eq. (flT)]) with 5=1 provided 

A 2 = l, 2A 3 = -A 5 [l + tanh 2 (^)] 1 A = 2tanh 2 (^)sech 2 (^) 

° L V M;J ' l + 2tanh 2 (/ l /3)-tanh 4 (^/3) V 7 

From Eq. (|3"6"j) it follows that A < A5/i 2 /2. Further, using the fact that tanh 2 (x) < 1, it is easily checked that the 
kink solution exists in this model only if < A < (2 — \/2)/2. 

In case only A3 and A5 are nonzero, the static Eq. {1} is equivalent to a highly nonlinear map which is even more 
complicated than the 18 parameter map of Quispel et al. [22l| . In fact, this is also true for the remaining three cases 
discussed below. 

The corresponding staggered solutions are easily found by changing the signs of terms with A 3 and (2 — A). In 
particular, it is easily shown that in the limit m = 1, the staggered sn solution reduces to the staggered kink solution 
Ijl9]) with 5 = (-1)" provided 

2A 3 = A 5 [l + tanh 2 ( ^)], A 5 h 2 = 2 -^— , A 2 = P + tnh 2 M tanh 2 ^) 

1 1 ;l 3 + tanh 2 (h/3) [2 + 3 tanh 2 (ft/3) - tanh 4 (/i/3)] V 7 



7 



< A < 4 . (38) 



One can show that this solution can only exist provided 

6 + V2 
2 

Another exact JEF solution is the nonstaggered dn solution as given by Eq. (fT8|) with 5=1 provided the following 
relations are satisfied: 

2 - A 

A 5 cs(2h/3, to) = -A 3 cs(h(3, to) , .„ TO = A 3 cs 2 (h/3, m) , (39) 

o 

A 3 ds(h(3, m)vs(hj3, to) + A 5 [ns(2h0, m)ds(2hj3, to) - cs 2 (2h0, to)] . (40) 



A 2 h 2 



Yet another exact JEF solution is the nonstaggered cn solution as given by Eq. (fTTj) with 5=1 provided the 
following relations are satisfied: 

A 5 ds{2hp,m) = -A 3 ds(hp,m), ^^^T = A 3 ds 2 (hp,m) , (41) 
= A 3 cs(h/3,m)ns(h/3,m) + A 5 [ns(2h@ , m)cs(2h/3 , m) - ds 2 (2h/3,m)} . (42) 

A z h z 

In the limit to = 1, both the cn and the dn solutions go over to the pulse solution Eq. ([20]) with 5=1, and the 
relations (|39p and (|4U|) as well as (|4"Tj) and (|4"2"|) take the simpler form 

A = -2[cosh(/i/3) - 1] < , A 5 <0, 2A 3 = -A 5 sech(/i/3) , |A|/i 2 = |A|j^| , 
A2 = (|A| + l)(|A| + 2)(|A|+4) ^ 



One can also obtain the corresponding staggered cn and dn solutions by changing the signs of ^3 and (2 — A) in 
the relations (f39|) and (j40| as well as |4T|) and (|42|) . In the limit to = 1 one can show that no staggered pulse solution 
Eq. (f2"0"l) with 5 = (—1)™ exists in case only A 3 and A 5 are nonzero. 

D. Case (v): only A%, A3, and A5 are nonzero 

Since only A%, A 3 , A$ are nonzero, hence Eq. (|2|) reduces to 

A 2 + A 3 + A b = X . (44) 
In this case, one has the sn solution Eq. (IT)]) with 5 = 1 provided the following relations are satisfied: 

(2 — A) 771 

A 5 ns(2/i/3, 777) = -A 3 ns(h/3, m) , ^ 2 ^ 2 = A 2 cs(ft,/3)ds(/i/3) - A 3 ns 2 (h/3, to) , (45) 

A2ns 2 (/7/3, to) — A3Cs(/7/3, m)ds(h/3, to) 



-A 5 [cs(2/i/3,m)ds(2/i/3,m) - ns 2 (2/i/3, to)] . (46) 



In the m = l case, for the kink solution Eq. (|T9|) with 5=1, the relations (|45j) to (|46|) take the simpler form 
A 2 = 1 , A7 2 A, 



. , , , , , _ [A - 2 tanh 2 (hp)} [1 + tanh 2 (h{3)} 



tanli (/i/3)[3-tanli (fr/3)] 

, [A + 2(A - 1) tanh 2 (ft/?) + (2 - A) tanh 4 (fe/3)] 2 [A - 2 tanh 2 (fe/3)] 

A2/1 = n 5 , A^h = t 5 . (47 

tanh 2 (/ty?)[3- tanh 2 (ft/3)] tanh 2 (ft/3) [3 - tanh 2 (ft/3)] 



Note that this solution is valid for any value of A including A = 2. 
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The corresponding staggered sn solution is easily obtained by changing the signs of terms with A 3 and (2 — A). In 
particular, it is easily shown that in the limit m = 1, the staggered sn solution reduces to the staggered kink solution 
(p]l with 5 = (-1)" provided 



h 2 A 3 
h 2 A 5 
h 2 A 2 



A[l + tanh 2 (/i/3)][A - 2(2 - tanh 2 (/i/3)] 
[Atanh 4 (/i/3) + (8 - A) tanh 2 (/i/3) + 2(A - 4)] ' 

2A[A-2(2-tanh 2 (/i/3)] 
[Atanh 4 (/i/3) + (8 - A) tanh 2 (/i/3) + 2(A - 4)] ' 
A[(A - 2) tanh 4 (/i/3) + 2(3 - A) tanh 2 (fe/3) + 4 - A] 

[Atanh 4 (/i/3) + (8 - A) tanh 2 (/i/3) + 2(A - 4)] 



a _ (A -4) 

^ - (2^ _A)' (48) 

Note that this solution exists over a vast range of A values including A = 2. 
For the dn solution Eq. (|18[) with 5=1 one has 

2 — A 

A 5 cs(2h/3, m) = -A 3 cs(h/3, m) , = A 3 cs 2 (h(3 , m) - A 2 ds(h(3, m)ns(h/3, to) , (49) 

o 

-A 2 cs 2 {h(3, to) + A 3 ds(hf3, m)ns(hp , to) + A§[ns(2hp , m)ds(2hf3 , to) — cs (2h(3, to)] . (50) 



A 2 h 2 

Parameters of the cn solution Eq. (fT7| with 5=1 satisfy 

(2 - A)to 



A 5 ds(2h(3,m) = -A 3 ds(hf3,m) , - . 2 " 2 / "" = A 3 ds 2 (fe/3, to) - A 2 ns(h/3,m)cs(h/3,m) , (51) 
= -^ 2 ds 2 (/i/3, to) + A 3 cs(/i/3, m)ns(h/3, to) + A 5 [ns(2/i/3, m)cs(2h(3, to) - ds 2 (2/i/3, to)] . (52) 



A 2 h 2 



In the limiting case of to = 1, from the relations (|4"9"| and (f50"|) as well as (4.20) and (4.21), one finds the relations 
for the pulse solution Eq. (|2"0|) with 5=1 

A = -2[cosh(/i/3) - 1] < , 2A3 = -A 5 sech(/i/3) , 

A 2 h 2 = ^[2coahW) - 1] - 2[cosh(^) - 1] ,A 2 = • (53) 

Thus the pulse solution exists only if A < 0. 

One can also obtain the corresponding staggered dn and cn solutions. In the limit to = 1, from these solutions we 
obtain the staggered pulse solution Eq. (f2"0")) with 5 = (—1)" satisfying 

A = 2[1 + cosh(fy3)] > 4, 2A 3 = A 5 sech(hf3)] , A 2 h 2 = -/i 2 A 3 [2cosh(/i/3) + 1] + 2[1 + cosh(/i/3)] , 
(A-2)(A-4) 

2[(A-3)ftM 3 -(A-2)] " 1 ' 

Thus the staggered pulse solution exists only if A > 4. 

E. Case (vi): only A3, A4, and A 5 are nonzero 

Since only A 3 , A4, A$ are nonzero, hence Eq. @ reduces to 

A 3 + A i + A 5 = \. (55) 
In this case, one has the sn solution Eq. (fl6|) with 5 = 1 satisfying 

(2 — Aim 

A 5 ns(2/i/3, to) = -A 3 ns(fy?, to) , 1 ; = (At - A 3 )ns 2 (/i/3, to) , (56) 
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2m 

„ = 2A4ns(hf3, m)ns(2h(3, m) — A 3 cs(/i/3, m)ds(h/3, m) 
-A 5 [cs{2h/3, m)ds(2/i/3, m) - ns 2 (2/i/3, m)} . (57) 



In the limit m = 1, we obtain from the sn solution the kink solution Eq. (| 19(1 with S = 1 with the parameters 
satisfying 



, 2 ,2, [l + tanh^/?)][(A-2)tanh^(fy?)+A] 
A = 1 , h A^ = g - 



2tanri (/i/3) 

^ = A + 2(A - 1) tanh 2 (/i/3) - (A - 2) tanh 4 (^) 2 ^ (A - 2) tanh 2 (fr/3) + A 

4 2tanh 2 (/i/3) ' 5 tanh 2 (fy3) ' 1 ' 

Note that this solution is valid for any A including A = 2 when conditions (|58|) take particularly simple form. 

One can also obtain the corresponding staggered sn solution. In the limit to — 1, we obtain from the staggered sn 
solution the staggered kink solution Eq. (fT9|) with 5 = (— 1)" provided 

A 2 = (3A-8)tanh 2 (fe/3) + 2(A~4) 
Atanh 2 (/i/?) 

4 h 2 = A[l + tanh 2 (/i/3)][4 - A + (2 - A) tanh 2 (fe/3)] 
3 ' 2[(8-3A)tanh 2 (/i/3) + 2(4-A)] ' 1 ' 

From the above equations it follows that no solution exists in case < A < 3. In particular, note that a solution does 
exist in case A < provided A 3 < 0. 

Parameters of the dn solution Eq. (fl~8|) with 5 — 1 can be found from 



2 - A 

A 5 cs(2h/3, m) = -A 3 cs(h{3, to) , = (A 3 - A 4 )ca 2 (hf3, to) , (61) 

2 

= —2AiCs(hp, m)cs(2h(3, m) + A 3 ds(/i/3, m)ns(hp, m ) + A b [ns(2/i/3, m)ds{2hf3, m) - cs 2 {2h/3, m)] . (62) 
On the other hand, the parameters of the cn solution Eq. (fT7)l with 5=1 satisfy 

A 5 ds{2h(3, to) = -A 3 ds(/i/3, m) , (2 ~ 2 ^ m = (A 3 - A 4 )ds 2 (/^/3, to) , (63) 

2TO s , „ , „ v „ . r „ x „ s ,2/ 



, „ = -2A 4 ds(ft,/3,TO)ds(2/i/3,TO) + A 3 cs(/i/3,TO)ns(/i/3,TO) + A 5 [ns(2/i/3,TO)cs(2/i/3,TO)-ds (2hf3,m)]. (64) 

In the limit to = 1, we obtain the parameters of the pulse solution Eq. (|20p with 5 = 1 

A = -2[cosh(fy3) - 1] < 0, 2A 3 = -A 5 sech(h(3) , A 4 h 2 = h 2 A 3 [2 cosh(/i/3) - 1] - 2[cosh(/i/3) - 1] , 
(|A|+2)(|A|+4) 

^ " 4(1 -W 3 ) ■ (65) 

Thus the pulse solution exists only if A < 0. 

One can also obtain the corresponding staggered dn and cn solutions. In the limit m = 1, from these solutions we 
obtain the staggered pulse solution Eq. (f2T)]) with 5 = (—1)™ satisfying 



A = 2[1 + cosh(/i/3)] > 4, 2A 3 = A 5 sech(/i/3) , A^ 2 = -/i 2 A 3 [2cosh(/i/3) + 1] + 2[1 + cosh(/i/?)] , 
^= '*- 2) '*- 4) . (66) 



Thus the staggered pulse solution exists only if A > 4. 
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F. Case (vii): nonzero A2, A3, A4, and As with Ai — Aa = 
In this case, Eq. ([2]) reduces to 

A 2 + A s + A4, + A 5 = A . (67) 
In this case, the sn solution Eq. (fT6|) with S = 1 is characterized by 

A s w(2hp, m) = -A 3 as(h(3, m) , (68) 

2m 

= 2Atns(/i/3, m)ns(2/i/3, m) — A3Cs(/i/3, m)ds(hj3, m) 



A 2 h 2 

+A 2 ns 2 (h(3) - A 5 [cs(2/i,5,m)ds(2/i/?,m) - ns 2 (2hf3,m)} , (69) 

(2 7 2 f 2 )W = (A 4 - A 3 )ns 2 (hf3, m) + A 2 cs(h/3)ds(h(3) . (70) 

In the limit m — 1, this reduces to the kink solution satisfying 

^ = 1, tanh 2 ^).-^^, Ai=l 2 + h2{2 ^ +A5) -A 3 -2A 5 . (71) 

Note that this solution is valid for any A including A = 2. 

One can also work out the corresponding staggered sn solution Eq. (fl"6|) with 5 = (—1)"- In the limit m — 1, we 
obtain the staggered kink solution Eq. (fTI)]) with S = (—1)™ satisfying 

A 2 = (4-A)[l + tanh 2 (^)] 

4A 3 /i 2 + (2A 2 /i 2 -A)[l + tanh 2 (/i/3)] ' 1 J 



2A 3 /i 2 [2(4 - A) + (8 - 3A) tanh 2 (/i/3)] 
[1 + tanh 2 (/i/3)] 



A 2 h 2 [2(4 - A) - Atanh 2 (/ l /3)] + ^ r / / , 2 , = A[4 - A + (2 - A) tanh 2 (/ l /3)] . (73) 



Another exact solution is the dn solution Eq. (fT8|) with 5 = 1 provided the following relations are satisfied: 

A 5 cs(2h/3, m) = -A 3 cs(h/3, m) , (74) 

2 

-2Atcs(/i/3, m)cs(2hf3, m) + A 3 ds(h(3, m)ns(/i/3, m) 



A 2 h 2 

-A 2 cs 2 (h(3) + A 5 [ns(2/i/3, m)ds{2h(3, m) - cs 2 {2hf3, m)] , (75) 

= (A 3 - A 4 )cs 2 (h(3, m) - A 2 ns(/i/3)ds(/i/3) . (76) 

Yet another exact solution is the cn solution Eq. (fT8"|) with S = 1 satisfying 

A 5 ds(2h/3, m) = -A 3 ds(h/3, m) , (77) 

2m 

.„, „ = — 2Atds(/i/3, m)ds(2hf3, m) + A 3 cs(hf3, m)ns(h(3, m) 
A z h 2 

-A 2 ds 2 {hf3) + A 5 [ns(2/i/3, m)cs(2h/3, m) - ds 2 (2/i/3, m)] , (78) 
(2 ^2fe 2 )W = (^3 ~ A 4 )ds 2 (/i/3, m) - A 2 ns(h(3)cs(hp) . (79) 
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In the limit m = 1, we obtain the pulse solution Eq. (|20[) with S = 1 

A = -2[cosh(>/?) - 1] < 0, 2A 3 = -A>sech(/i/3) , (A* + A 2 )h 2 = h 2 A 3 [2 cosh(/i/3) - 1] - 2[cosh(fy3) - 1] , 
2 _ (|A|+2)(|A|+4) " 

2(2 - 2h 2 A 3 - A 2 h 2 ) ' 1 ' 

Thus the pulse solution exists only if A < 0. 

One can also work out the corresponding staggered dn and cn solutions. In the limit m = 1, we obtain the staggered 
pulse solution Eq. $2$ with S = (-1)" with 

A = 2[1 + cosh(/i/3)] > 4, 2A 3 = A 5 sech(h/3) , (Ai + A 2 )/i 2 = -h 2 A 3 [2cosh(h/3) + 1] + 2[1 + cosh(ft/3)] , 
(A-2)(A-4) 

2(2 - 2^ 2 A 3 - A 2 /i 2 ) ' 1 ; 

Thus the staggered pulse solution exists only if A > 4. 

Summarizing, it is worth noting that in all the five cases discussed above (which is also true in the remaining two 
cases (i) and (ii)), for the kink solution, the amplitude A is always equal to 1. Further, in all the seven cases, while 
the pulse solution exists only if A < 0, the staggered pulse solution exists only if A > 4. 

It is worth pointing out here that the trigonometric solutions for these models do not follow from above JEF 
solutions in the limit m = 0, since both left and right hand sides of the identities for the Jacobi elliptic functions sn 
and cn, which one has used in deriving these solutions, vanish identically in this limit (23j . However, trigonometric 
solutions can be derived for these models directly, as we show below. 



G. Trigonometric solutions 

We shall now show that unlike the JEF and the hyperbolic solutions, the static TI trigonometric solutions with an 
arbitrary shift along the chain xq exist even when all the six parameters A are nonzero. Ideally, the trigonometric 
solutions should be the m — ► limit of the JEF solutions. However, the relevant JEF identities we have used vanish 
in this limit and thus the trigonometric solutions must be derived separately. 

We look for a solution of Eq. {I) of the form 

4> n = SA sm[hf3(n + x )] , (82) 

where S = 1 or (— 1)™ depending on whether it is a nonstaggered or a staggered solution, and find that it exists under 
the following two conditions 

± (A - 2) + 2 cos(hj3) = h 2 A 2 sm 2 (hj3)[±(A 3 - At) + (3A 6 - As) cos(h(3)] , (83) 



± [A + At + As cos(2h/3)] + (A 2 + As) cos(/j/3) + As cos(/i/3)[4cos 2 (/i/3) - 3] = , (84) 

where the upper (lower) sign corresponds to nonstaggered (staggered) sine solution. As mentioned before, if 4> n (t) 
is a solution to Eq. |T]), then the staggered solution (— l) n (f) n (t) is also a solution to Eq. |T]) with the coefficients of 
A\, A3 and A4 having opposite signs and with 2 — A replaced by A — 2, where A = h 2 X. To make the presentation 
compact, in most cases we shall therefore only give results for the nonstaggered solutions. We shall see that these 
solutions exist in all the models discussed in this paper except in the case of Model 1 (only A\ nonzero) and Model 4 
(only A4 nonzero). 

Using the well known addition theorem for sine, it is easily shown that the sine solution follows from the two-point 
quadratic map 

4>l+i + & - ^K+iK cos{hf3) ~ A 2 sin 2 (h0) = . (85) 

We now discuss these solutions in the several special cases. 
Model 1: Only A ^ 0, {h 2 A x = A) 

From Eqs. (f8"3"|) and (|5^|) it follows that in this case there is neither a staggered nor a nonstaggered solution of the 
form of Eq. l(82]). 

Model 6: Only Ai ^ , (h 2 A 4 = A) 

From Eqs. (|8"3"|) and (IM)) it follows that in this case also there is neither a staggered nor a nonstaggered solution of 
the form of Eq. (|g2"]). 
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Model 3: Only A 2 ^ , (h 2 A 2 = A) 

In this case, both staggered and nonstaggered solutions exist provided 

A = 2, = ±7r/2. (86) 

and in this case A is arbitrary. 

Only A 4 , A 2 ^0, [h 2 (A 4 + A 2 ) = A] 

In this case, the nonstaggered solution with 5 = 1 given by Eq. (l82|) exists provided 

A A + A 2 cos(fy3) = , 2A 2 Acos(ft / 3)cos 2 (h/3/2) = A - 2[1 - cos(/i/3)] . (87) 

On the other hand, the staggered solution with S = (—1)" exists in the same model (with only A 2 and At nonzero) 
provided 

Ai + A 2 cos(ftj3) = , 2A 2 Acos(/i/3) sin 2 (ft/3/2) = A - 2[1 + cos(/i/?)] . (88) 

Only A 4 , A 3 , A 5 ^0, [h 2 (A 4 + A 3 + A 5 ) = A] 

In this case the nonstaggered solution exists provided 

A± + A 5 cos(/i/3) + A 3 cos(2hp) = , 2A 2 h 2 A 3 cos 2 {hf3) sin 2 (ft/3) = A — 2[1 — cos(ft/3)] . (89) 

Only A 2 , A 3 , A 5 + 0, [ft 2 (A 2 + A 3 + A 5 ) = A] 

In this case the nonstaggered solution exists provided 

(A 2 + A 5 ) cos(hp) + A 3 cos(2ft/3) = , A 2 h 2 [A 3 - A 5 cos(h/3)} sin 2 (ft/?) = A — 2[1 - cos(ftjS)] . (90) 

Only A 3 ,A 5 ^Q [h 2 (A 3 + A 5 ) = A] 

In this case the nonstaggered solution exists provided 

A 5 cos(hP) + A 3 cos(2ft/3) = , 4A 2 A cos 3 (hp) cos 2 (^) = (A - 2[1 - cos(fc/3)])[l + 2 cos(fy3)] . (91) 

Case of Ai = A 6 = 0, A 2 , A 3 , A 4 , A 5 ^ , [/i 2 (A 4 + A 2 + A 3 + A 5 ) = A] 
In this case the nonstaggered solution exists provided 

A A + {A 2 + A 5 ) cos(/i/3) + A 3 cos(2ft/3) = , A 2 h 2 [A 3 - A4- A 5 cos(hp)} sin 2 (/i/3) = A — 2[1 — cos(/i/3)] . (92) 

Model 2 with 5 = 0,7 = 1/4: Only A 5 = A 6 7^ , 2ft 2 A 5 = A 
In this case the nonstaggered solution exists provided 

cos 2 (hp) = - , A 2 A = 4 ± 2 V2(A - 2) . (93) 

Thus in this case one has a solution with period 8. 

Model 2 with 7 = 0,5 = 1/4: Only A\ = A 3 /2 = A 4 = A/2 are nonzero, 2h 2 A 3 = A 
In this case the nonstaggered solution exists provided 

sin 2 (hp) = 1, A 2 A = 4(A-2). (94) 

Thus in this case one has a solution with period 4. 

Model 2: Case of A 3 = 2A X = 2A A , A 6 = A 5 ,A 2 all nonzero , [h 2 (2A 3 + 2A 5 + A 2 ) = A] 
In this case the nonstaggered solution exists provided 

cos(hp)[A 2 + 2A 5 cos(2hp) + 2A 3 cos(ft/3)] = , A 2 h 2 [A 3 + AA 5 cos(ft/3)] sin 2 (ft/3) = 2 [A - 2(1 - cos(ft/3))] . (95) 

Thus in this case one has solution with period 4. 

Model 10: Hamiltonian model given by Eq. (|15|) with A% = Act\\ , A 2 = 6a 2 A , A 3 — 4a 3 X , A$ — 2a 2 X , A4 = 
A 5 =0 [A(a t + 2a 2 + a 3 ) = 1] 

In this case the nonstaggered solution exists provided 

2A 2 A[2a 2 + 3a 2 cos(ft/3)] sin 2 (hp) = [A - 2(1 - cos(ft/3))] , 

a 2 cos 3 (hp) + a 3 cos 2 (ft/3) = a 2 + a 3 — i . (96) 
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Speight and Ward Model 7: ai = a 2 = §013 = 

Being a special case of the model 10, in this case Eq. reduces to 

4cos 3 (/i/3) +6cos 2 (/i/3) -1 = 0, A 2 A[1 + cos(h(3)\ sm 2 (h/3) = 3[A-2(1 - cos(h/3))} . (97) 

The roots of the first equation in Eq. (|97p are h(3 — 2ir/3 and h(3 « 1.196. The first root corresponds to the three- 
periodic TI solution with an amplitude that can be found from the second equation in Eq. (|97|) as A 2 — 8(A — 3)/ A, 
so that the solution is defined either when A > 3 or when A < 0. The second root corresponds to the TI solution with 
the period approximately 5.253. 

Instead, the staggered solution exists in this model provided 

4cos 3 (/i/3) -6cos 2 (fy3) + 1 = 0, A 2 A[1 - cos(/i/3)] sin 2 (/i/3) = 3 [A - 2(1 + cos(/i/3))] . (98) 

The roots of the first equation in Eq. |98|) are h(3 = ir/3 and hf3 « 1.944. The first root corresponds to the six-periodic 
TI solution with the amplitude which is same as for the three-periodic nonstaggered solution, i.e. A 2 = 8(A — 3)/A, 
so that this solution is also defined either when A > 3 or when A < 0. The second root corresponds to the TI solution 
with the period approximately 3.023. 

It is worth pointing out that in the case of Model 3 (only A2 nonzero), both the staggered and the nonstaggered 
solutions are of period 4, i.e. h(3 = ir/4. The same is also true of Model 2 in case S = 1/4,7 = 0, i.e. Ai — A4 = 
A3/2 = A/4. On the other hand, in the case of Model 2 with 7 = 1/4, S = 0, i.e. A^ — Aq ~ A/2, both the staggered 
and the nonstaggered solutions are of period 8. Finally, in the case of Model 7, i.e. Speight and Ward model with 
Ai = 2A/9, A 2 = A 3 = A/3, A4 = A 5 = Q,A 6 = A/9, one of the nonstaggered solution is three-periodic while in the 
staggered case, one of the solution is of period 6. 



IV. SHORT-PERIOD SOLUTIONS 



We shall now show that apart from the JEF, hyperbolic and trigonometric, there are also several short period 
and even aperiodic solutions of Eq. (JTJ) . In order to obtain these solutions, it is useful to look at the symmetries of 
Eq. ([1]). In particular, notice that Eq. ([I]) is invariant under </> n _i — > <fi n +i and <fi n +i — > 4> n -i- Further, Eq. ^ 
is also invariant under (0 n _i, <f> n +i) —* (— </>n-ij ~4 > ni —4>n+i)- A consequence of these two symmetries is that if 
((f>n-i, <t>m <t>n+i) is a solution to Eq. (P) under certain constraints, then {—<j> 

l), (<t>n+l,<t>n,<t>n-l) and 

(— <t> n +x, ~4> n , —<t>n-i) are also solutions of Eq. (QJ provided the same constraints are satisfied. 

We list below several exact solutions to Eq. ([I]) satisfying the constraint ((2]). We shall only write down the exact 
solutions to the general model. IN Appendix we spell ot the short-period solutions that are admitted by many of the 
models discuused in this paper. 

While obtaining the periodic solutions the following results have been used which have been derived by using Eqs. 
© and ©• 

1. If 4>n-i = 4>n = 4>n+i = a then it follows that 

a 2 = 1 . (99) 

2. If 4> n -\ — cf> n = —(f>n+i = o, then it follows that 

A-2 = h 2 a 2 (A 1 +A 3 -A 4 ). (100) 

3. If 4> n -\ = —(fin = <fin+i = o> then it follows that 

A - 4 = /iV(Ai -A 2 + A 3 + A 4 - A 5 - A 6 ) . (101) 

4. If (fi n -i — <fi n — a, and <fi n +i = then it follows that 

2(A - 1) = h 2 a 2 {2Ai + A 2 + A 3 + A 6 ) . (102) 

5. If (fin-i = 4>n+i = a, and <fi n = then it follows that 

2 = h 2 a 2 {A 5 + A e ). (103) 
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6. If 4> n -i — —<p n = a, and 4>n+i — then it follows that 

2{A-3) = h 2 a 2 (2Ai-A 2 + A 3 -A e ). (104) 

7. If 4> n -i = a, and 4> n = 4>„.+i = then it follows that 

2 = h 2 a 2 A 6 . (105) 

8. If <j) n — a, and </> n _i = 4> n +i = then it follows that 

A - 2 = h 2 a 2 A 1 . (106) 

We now discuss the various exact solutions. 

(i) Solution with period 2: <p = (■■■, a , —a, ■■■) 

This is an exact solution to Eq. ([T]) satisfying the constraint {2} provided 

A- 4 = h 2 a 2 \A x -A 2 +A 3 + A 4 -A 5 - A 6 ] . (107) 

In case A = 4, then one has a solution with a being an arbitrary real number provided h 2 (A\ + A3 + A 4 ) = 
h 2 (A 2 + A 5 + A 6 ) = 2. 

(ii) second solution with period 2: <f> — (..., a, 0, ...) 

This is an exact solution to Eq. ([1]) satisfying the constraint ([2]) provided 

A-2 = h 2 a 2 A 1 , 2 = (A 5 + A 6 )h 2 a 2 . (108) 

One can in fact generalize this solution and show that even 
4> =(..., a, 0, (a, p times), — a, 0...) 

is an exact solution with period 2p + 2 (p > 2) provided Eq. (|108|) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 2p + 2) and 
with "—a, 0" added randomly after any "0" is an exact solution to Eq. (fT]) provided Eq. I|108p is satisfied. 

(iii) third solution with period 2: (f> = (..., a, ~, a, -, ...) 

This is an exact solution to Eq. (fT]) satisfying the constraint ^ provided 

A 1= A 5 + A 6 =Q, h 2 A 2 =A-2, h 2 (A 4 + A 3 ) = 2, (109) 

and in this case a is an arbitrary real number. 

(iv) Solution with Period 3: cf) = (..., a, —a, 0, ...) 

It is easily shown that this is an exact solution to Eq. |T|) satisfying the constraint ([2]) provided 

a 2 h 2 [2A 1 + A 3 -A 2 -A 6 ] = 2(A-3). (110) 

In case A = 3, then one has a solution with a being arbitrary real number provided 2Ai + A3 = A 2 + Ag. 

(v) second solution with Period 3: <f> = (..., a, a, —a, ...) 

It is easily shown that this is an exact solution to Eq. {T]) satisfying the constraint ([2]) provided 

a 2 h 2 {A 1 + A 3 - A 4 ) = A -2, a 2 h 2 (A 1 + A3 + A 4 - A 2 - A 5 - A e ) = A - 4 . (Ill) 

One can in fact generalize this solution and show that even 
<j} = (..., a, —a, (a, —a p times), —a,...) 

is an exact solution with period 2p+l (p > 1) provided Eq. (|111[) is satisfied. 

One can also show that an aperiodic solution with any number of "a" and "—a" kept at random but with the 
constraint that at most two "a" or two "—a" are always together, is an exact solution to Eq. (fT]) provided Eq. pil[) 
is satisfied. 

(vi) third solution with Period 3: (j> = (..., a, a, 0, ...) 

It is easily shown that this is an exact solution to Eq. |T]) satisfying the constraint ([2]) provided 

a 2 h 2 (A 2 + 2A 1 +A 3 + A 6 ) = 2(A-l), h 2 a 2 (A 5 + A 6 ) - 2 . (112) 

One can in fact generalize this solution and show that even 

4> = (..., a, a, 0, —a, —a, 0, (a, a, 0, —a, —a, p times), a, a, 0, ...) 

is an exact solution with period 6p + 3 (p > 1) provided Eq. (|112p is satisfied. 
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One can also show that an aperiodic solution can be constructed from the above periodic solution with period 6p + 3 
with "a, a,0" or "-a, -a,0" added at random between "0" and "-a" or "0" and "a". 

(vii) Solution with Period 4: <fi — (..., a, a, a, 0, ...) 

It is easily shown that this is an exact solution to Eq. (JTJ) satisfying the constraint (J5]) provided 

a 2 = l, h 2 {A 2 + 2A 1 + A 3 +A 6 ) = 2(A- 1) , /i 2 (A 5 + A 6 ) = 2 . (113) 

From here it follows that such a solution is valid provided A 2 + A3 + 2A4 + A$ =0. 
One can in fact generalize this solution and show that even 
4> =(..., a, a, 0, (a, a, p times), a,...) 

is an exact solution with period 3p+ 1 (p > 1) provided Eq. (|113|) is satisfied. 

One can also show that an aperiodic solution with any number of "a" and "0" but with the constraint that at least 
two "a" are always together and no two "0" are either nearest or next-to- nearest neighbours is also an exact solution 
provided Eq. (| 1 1 3[) is satisfied. 

(viii) second Solution with Period 4: <fr = (..., a, a, a, — a, ...) 

It is easily shown that this is an exact solution to Eq. (JXJ) satisfying the constraint @ provided 

a 2 = l, h 2 (A 1 +A 3 -A 4 )=A-2, h 2 {A x + A 3 + A 4 ~ A 2 - A 5 ~ A 6 ) = A - 4 . (114) 

In view of the constraint ^ this implies that such a solution is valid only if A4 = 0. 
One can in fact generalize this solution and show that even 
4> =(..., a, a, a, (a p times), —a,...) 

is an exact solution with period p + 1 (p > 3) provided Eq. (|114[) is satisfied. 

One can also show that an aperiodic solution with any number of "a" and "—a" but with the constraint that at 
least once three or more "a" or "—a" are together. 

(ix) third solution with Period 4: <p — (■••■ a , b, —a, —b, ...), where a 2 7^ 6 2 . 

It is easily shown that this is an exact solution to Eq. (JTJ) satisfying the constraint @ provided 

(a 2 + b 2 )h 2 Ai = A — 2, A3 = A a + A 4 , a 2 7^ 6 2 , Ai^O. (115) 

Thus one has a one parameter family of solutions. It is easily seen that such a solution will exist in model 2 (in case 
5 is arbitrary but nonzero, while 7 is arbitrary), model 8 (in case 7 = 0, a = 1/4) and Hamiltonian model 10 (in case 
ai = a 3 ). 

In the special case of A 3 = A4, A\ = and A = 2, one, in fact, has a two parameter family of solutions in the sense 
that now both a and b are arbitrary real numbers. 

(x) fourth solution with Period 4: = (..., a, a, — a, — a, ...) 

It is easily shown that this is an exact solution to Eq. |T|) satisfying the constraint @ provided 

a 2 h 2 [A 1 + A 3 - Ai] = A- 2. (116) 

In case A = 2, then one has a solution with a being an arbitrary real number provided A\ + A 3 — A4. 

(xi) fifth solution with Period 4: = (..., a, 0, — a, 0, ...) 

It is easily shown that this is an exact solution to Eq. (JTJ) satisfying the constraint ^ provided 

a 2 h 2 A 1 = A- 2. (117) 

In case A = 2, then one has a solution with a being an arbitrary real number provided A\ = 0. 

(xii) sixth solution with Period 4: (j> = (..., a, 0, —a, a, ...) 

It is easily shown that this is an exact solution to Eq. (JTJ) satisfying the constraint @ provided 

a 2 h 2 [A 1 +A 3 - A A ] = A- 2, a 2 h 2 [2A l + A 3 + A 2 + A 6 ] = 2(A- 1) , a 2 h 2 [2A l + A 3 - A 2 - A 6 ] = 2(A - 3) . (118) 

From here it follows that such a solution is valid provided A 3 = 2A4. 
One can in fact generalize this solution and show that even 
4> =(..., a, 0, —a, (a, 0, —a p times), a,...) 

is an exact solution with period 3p+ 1 (p > 1) provided Eq. (|118[) is satisfied. 

One can also show that an aperiodic solution with any number of "a" and "—a" but with the constraint that at most 
two "a" or two "—a" are always together and further "0" are added at random between "—a" and "a" or between "a" 
and "—a" but with the proviso that no two "0" are ever nearest or next-to-nearest neighbours, is an exact solution 
to Eq. (|T|) satisfying the constraint ^ provided Eq. (| 1 1 8|) is satisfied. 

(xiii) seventh solution with Period 4: <fi = (..., a, — a, 0, a, ...) 
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It is easily shown that this is an exact solution to Eq. |T|) satisfying the constraint ([2]) provided 

a 2 h 2 [A 1 +A 3 + A 4 -A 2 -A 5 -A 6 }=A-4, cfh 2 ^ + A 3 - A 2 - A e ] = 2(A - 3) , a 2 h 2 [A e + A 5 ] = 2 . (119) 

From here it follows that such a solution is valid provided A 2 + Ag = A 3 + 2A 4 . 
One can in fact generalize this solution and show that even 
<j> = (..., a, —a, 0, (a, —a, 0, p times), —a,...) 

is an exact solution with period 3p + 1 (p > 1) provided Eq. (|119[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 3p+ 1) by 
adding at random "a, —a" between any and "a" is an exact solution to Eq. {T]) satisfying the constraint ((2]) provided 
Eq. ([Tig]) is satisfied. 

(xiv) Solution with Period 5: cf) — (..., a, a, a, —a, —a, ...) 

It is easily shown that this is an exact solution to Eq. |T]) satisfying the constraint ([2| provided 

a 2 = l, h 2 \A x + A 3 - Aa] = A - 2. (120) 

One can in fact generalize this solution and show that even 
cj> = (•■•, a, a, a, (a p times), — a, — a, (— a q times), ...) 

is an exact solution with period p + q (p > 3, q > 2 or p > 2, q > 3) provided Eq. (|120[) is satisfied. 

One can also show that an aperiodic solution constructed from any numbers of "a" and "—a" with the constraint 
that two or more of "a" as well as two or more of "—a" are always together is an exact solution to Eq. (p} satisfying 
the constraint ^ provided Eq. (|120|) is satisfied. 

(xv) second solution with period 5: <fi = (..., 0, a, 0, a, a, ...) 

It is easily shown that this is an exact solution to Eq. |T]) satisfying the constraint ([2| provided 

h 2 a 2 [2A 1 + A 3 + A 2 + A 6 ] = 2(A - 1) , h 2 a 2 (A 6 + A 5 ) = 2 , h 2 a 2 A 1 = A - 2 . (121) 

From here it follows that such a solution is valid provided A 2 + A 3 = A$ . 
One can in fact generalize this solution and show that even 
(f> = (...0,a, (0, a p times), a,...) 

is an exact solution with period 2p + 1 (p > 2) provided Eq. (I121|l is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 2p+ 1) by 
adding at random "—a, — a,0" after any "0" is an exact solution to Eq. |T]) satisfying the constraint @ provided Eq. 
(fT2Tj) is satisfied. 

(xvi) third solution with period 5:0= (..., a, —a, a, —a, 0, ...) 

This is an exact solution to Eq. ([T]) satisfying the constraint ([2]) provided 

A-4 = h 2 a 2 (A i + A 1 +A 3 -A 2 -A 5 -A 6 ), 2(A - 3) = h 2 a 2 (2A 1 + A 3 - A 2 - A 6 ) . (122) 

One can in fact generalize this solution and show that even 
4> = (...a, — a, (a, —a p times), 0,...) 

is an exact solution with period 2p + 1 (p > 2) provided Eq. (1122|) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 2p + 1) 
by adding at random "0" between "—a" and "a" or "a" and "—a", is an exact solution to Eq. ([T]) satisfying the 
constraint ([2]) provided Eq. (| 1 22|) is satisfied. 

(xvii) fourth solution with period 5: <fr = (..., a, 0, a, 0, — a, ...) 

This is an exact solution to Eq. ([T]) satisfying the constraint ([2]) provided 

A- 2 = h 2 a 2 A 1: 2 = h 2 a 2 (A 5 + A 6 ) , 2(A - 3) = h 2 a 2 (2A 1 + A 3 - A 2 - A 6 ) . (123) 

From here it follows that such a solution is valid provided A 2 = A 3 + A$ . 
One can in fact generalize this solution and show that even 
4> =(..., a, 0, (a, Op times), —a,...) 

is an exact solution with period 2p + 1 (p > 2) provided Eq. (I123|) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 2p + 1) by 
adding randomly "—a" between "0" and "a" with a constraint that no two "—a" can be either nearest or next-to- 
nearest neighbors, is an exact solution to Eq. ([I]) satisfying the constraint ([2]) provided Eq. (|123jl is satisfied. 

(xviii) fifth solution with period 5: <fi = (..., a, a, — a, — a, 0, ...) 

This is an exact solution to Eq. ([TJ) satisfying the constraint provided 



A - 2 = h 2 a 2 {A 1 +A 3 - A 4 ) , 2(A - 1) = h 2 a 2 {2A 1 +A 3 + A 2 + A 6 ) 



(124) 
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One can in fact generalize this solution and show that even 
(f> = (..., a, a, —a, —a, (a, a, —a, —a p times), 0,...) 

is an exact solution with period Ap+1 (p > 1) provided Eq. (|124|1 is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period Ap + 1) 
by randomly adding "0" between "—a" and "a" or "a" and "—a" with a constraint that no two "0" can be either 
nearest or next-to-nearest neighbors, is an exact solution to Eq. (p} satisfying the constraint ((2]) provided Eq. (jl24[) 
is satisfied. 

(xix) sixth solution with period 5:0 = (..., a, a, 0, —a, 0, ...) 

This is an exact solution to Eq. ([I]) satisfying the constraint ([2]) provided 

A- 2 = h 2 a 2 Ai , 2(A-l) = h 2 a 2 (2A 1 + A 3 + A 2 + A 6 ). (125) 

One can in fact generalize this solution and show that even 
<fi = (•■•, a, a, —a, —a, (a, a, —a, —a p times), 0,...) 

is an exact solution with period Ap+1 (p > 1) provided Eq. (I125|) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period Ap + 1) 
by randomly adding "0" between "—a" and "a" or "a" and "—a" with a constraint that no two "0" can be either 
nearest or next-to-nearest neighbors, is an exact solution to Eq. |T]) satisfying the constraint ((2]) provided Eq. (|125p 
is satisfied. 

(xx) seventh solution with period 5: (/) = (..., a, —a, 0, —a, a, ...) 

This is an exact solution to Eq. ([I]) satisfying the constraint ([2]) provided 

A- 2 = h 2 a 2 (A 1 + A 3 - A 4 ) , 2 = h 2 a 2 (A 5 + A 6 ) , 2(A - 3) = h 2 a 2 (2A 1 + A 3 - A 2 - A 6 ) . (126) 

From here it follows that such a solution is valid provided A 2 + A 3 = 2A 4 + A 5 . 
One can in fact generalize this solution and show that even 
(f> = (..., a, —a, 0, (a, —a, p times), —a, a,...) 

is an exact solution with period 3p+ 2 (p > 1) provided Eq. (I126P is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 3p + 2) by 
randomly adding "—a, a" between "0" and "a" is an exact solution to Eq. (p} satisfying the constraint |(2J provided 
Eq. (|T25)| is satisfied. 

(xxi) eighth solution with period 5: <fr = (..., a, 0, — a, a, a, ...) 

This is an exact solution to Eq. (fT]) satisfying the constraint ^ provided 

a 2 = l, A-2 = h 2 (A 1 +A 3 -A 4 ), 2(A- 3) = h 2 (2A 1 + A 3 - A 2 - A 6 ) , 2(A- 1) - /i 2 ^ + A 3 + A 2 + A 6 ) . (127) 

In view of the constraint such a solution is thus valid only if A 3 — 2A 4 = — A 5 . 
One can in fact generalize this solution and show that even 
<p =(■■■, a, 0, —a, (a, 0, —a p times), a, a,...) 

is an exact solution with period 3p+ 2 (p > 1) provided Eq. (|12T|) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 3p + 2) by 
randomly adding "a, a" between "—a" and "a" is an exact solution to Eq. ^ satisfying the constraint |(2} provided 
Eq. ([12"?)) is satisfied. 

(xxii) ninth solution with period 5: 4> = (..., a, a, 0, a, — a, ...) 

This is an exact solution to Eq. dJ) satisfying the constraint ^ provided 

2 = h 2 a 2 {A 5 + A 6 ), A-2 = h 2 a 2 {A 1 +A 3 -A 4 ), 2(A - 3) = h 2 a 2 {2A 1 + A 3 - A 2 - A 6 ) , 

2(A-l) = h 2 a 2 {2A 1 +A 3 + A 2 + A 6 ), A - 4 = h 2 a 2 {A 1 + A 3 + A 4 - A 5 - A 6 - A 2 ) . (128) 

From here it follows that such a solution is valid provided A 2 = A 3 = 2A 4 = Ag. 
One can in fact generalize this solution and show that even 
4> =(..., a, a, 0, (a, a, p times), a, —a,...) 

is an exact solution with period 3p+ 2 (p > 1) provided Eq. (|128[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 3p + 2) by 
randomly adding "—a" between two "a" is an exact solution to Eq. |T]) satisfying the constraint (J2J) provided Eq. 
(fT2"8")) is satisfied. 

(xxiii) Solution with period 6: <f> = (..., a, a, 0, —a, —a, 0, ...) 

This is an exact solution to Eq. (fl]) satisfying the constraint provided 



2(A - 1) = h 2 a 2 (2A 1 +A 3 + A 2 + A 6 ) 



(129) 
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In case A = 1, then one has a solution with a being an arbitrary real number provided 2Ai + A 3 + A 2 + Aq = 0. 

(xxiv) second solution with period 6: <\> = (..., a, —a, 0, —a, a, 0, ...) 
This is an exact solution to Eq. ([I]) satisfying the constraint J5J provided 

2(A- 3) = h 2 a 2 {2A 1 + A 3 - A 2 -A3), 2 = h 2 a 2 (A 5 + A 6 ) . (130) 

One can in fact generalize this solution and show that even 
4> =(..., a, —a, 0, (a, —a, p times), — a, a, 0, ...) 

is an exact solution with period 3p + 3 (p > 1) provided Eq. (|130[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 3p + 3) by 
randomly adding "—a, a, 0" between "0" and "a" is an exact solution to Eq. {TJ satisfying the constraint (J2) provided 
Eq. ([HQ]) is satisfied. 

(xxv) third solution with period 6: </>=(..., a, a, a, 0, —a, —a, ...) 
This is an exact solution to Eq. ([T]) satisfying the constraint (J5J) provided 

a 2 = l, 2(A- 1) = h 2 {2A 1 +A 3 + A 2 + A e ) , A - 2 = h 2 (Ai + A 3 - At) . (131) 

In view of the constraint ([2]) this implies that such a solution is valid only if A 3 = —A$. 
One can in fact generalize this solution and show that even 
(f> = (..., a, a, a, (a p times), 0, — a,— a,...) 

is an exact solution with period p + 3 (p > 3) provided Eq. (|131[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period p + 3) 
by randomly adding "0, —a, —a" between any two "a" is an exact solution to Eq. ([T]) satisfying the constraint ([2]) 
provided Eq. (|131[) is satisfied. 

(xxvi) fourth solution with period 6: = (..., a, a, 0, —a, 0, a, ...) 
This is an exact solution to Eq. (TT]) satisfying the constraint ^ provided 

a 2 = l, 2(A- 1) = h 2 (2A 1 + A 3 + A 2 + A3) , A-2 = h 2 A 1 . (132) 

In view of the constraint ^ this implies that such a solution is valid only if A^ = — A5. 
One can in fact generalize this solution and show that even 
4> =(..., a, a, 0, —a, 0, (a, a, 0, — a, p times), a,...) 

is an exact solution with period 5p+ 1 (p > 1) provided Eq. (|132p is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 5p + 1) 
by randomly adding as many "a" as one wants between any two "a" is an exact solution to Eq. ([1]) satisfying the 
constraint @ provided Eq. (|132[) is satisfied. 

(xxvii) fifth solution with period 6: <j) = (..., a, —a, a, a, —a, 0, ...) 
This is an exact solution to Eq. ([T]) satisfying the constraint provided 

2(A-3) = h 2 a 2 {2A 1 +A 3 -A 2 -A 6 ) , A-2 = h 2 a 2 {A 1 +A 3 - A 4 ) , A-4 = /iV^i+Aj + At- A 2 - A3- A3) • (133) 

From here it follows that such a solution is valid provided A 3 = A5. 
One can in fact generalize this solution and show that even 
4> =(..., a, —a, a, a, —a, (a, —a, a, a, —a p times), 0,...) 

is an exact solution with period 5p+ 1 (p > 1) provided Eq. (|133|) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 5p + 1) by 
randomly adding "—a" between any two "a" is an exact solution to Eq. ([T]) satisfying the constraint ^ provided Eq. 
(I133|) is satisfied. 

(xxviii) sixth solution with period 6: <fi = (..., a, a, 0, a, — a, 0, ...) 

This is an exact solution to Eq. (TT|) satisfying the constraint ^ provided 

2(A-l) = / l V(2A + A5 + A 2 + A ; ), 2(A-3) ^ h 2 a 2 (2A 1 + A 3 - A 2 - A3), 2 = h 2 a 2 {A 5 + A 6 ) . (134) 

From here it follows that such a solution is valid provided A 2 = A5. 
One can in fact generalize this solution and show that even 
4> =(..., a, a, 0, (a, a, p times), a, — a,0, ...) 

is an exact solution with period 3p + 3 (p > 1) provided Eq. (|134[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 3p + 3) by 
randomly adding "a, — a,0" between "0" and "a" is an exact solution to Eq. (TTJ) satisfying the constraint @ provided 
Eq. (TIM)) is satisfied. 



19 



(xxix) seventh solution with period 6: = (..., a, a, — a, 0, — a, a, ...) 
This is an exact solution to Eq. (Q]) satisfying the constraint © provided 

a 2 = l, (A-2) = h 2 (A 1 +A 3 ~A 4 ), 2(A - 3) = h 2 (2A 1 + A 3 - A 2 - A 6 ) , 2 = h 2 (A 5 +A 6 ). (135) 

In view of the constraint © this implies that such a solution is valid only if A 2 — ~2A 4 , A^ ~ 2A 2 + A 3 . 
One can in fact generalize this solution and show that even 
<j> = (..., a, a, —a, 0, —a, (a, a, —a, 0, —a p times), a,...) 

is an exact solution with period 5p + 1 (p > 1) provided Eq. (|135p is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 5p+ 1) by 
randomly adding as many "a" as one wants between two "a" , is an exact solution to Eq. (fTJ) satisfying the constraint 
© provided Eq. (fH5)) is satisfied. 

(xxx) eighth solution with period 6: <f> = (..., a, 0, a, 0, a, —a, ...) 
This is an exact solution to Eq. (JTJ) satisfying the constraint © provided 

A - 2 = h 2 a 2 A 1 , (A - 4) = h 2 a 2 (A 1 + A 3 + A 4 - A 5 - A e - A 2 ) , 

2(A-3) = h 2 a 2 (2A 1 +A 3 -A 2 -A 6 ), 2 = h 2 a 2 {A 5 + A 6 ) . (136) 

From here it follows that such a solution is valid provided A4 = A$ , A 2 = A 3 + A 4 . 
One can in fact generalize this solution and show that even 
<fi = (..., a, 0, a, 0, (a, 0, p times), a, —a,...) 

is an exact solution with period 2p + 2 (p > 2) provided Eq. (|136p is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 2p + 2) by 
randomly adding "—a" between "0" and "a", is an exact solution to Eq. {1} satisfying the constraint ([2|) provided 
Eq. ([Hp]) is satisfied. 

(xxxi) ninth solution with period 6: <p = (..., a, 0, a, 0, a, a, ...) 

This is an exact solution to Eq. satisfying the constraint ((2j) provided 

A - 2 = h 2 A l , (A - 2) = h 2 {A l + A 3 - A A ) , a 2 = 1 , 2 = h 2 {A b + A 6 ) . (137) 

In view of the constraint © this implies that such a solution is valid only if A 2 = —2A 3 = —2A4. 
One can in fact generalize this solution and show that even 
cj) = (..., a, 0, a, 0, (a, 0, p times), a, a,...) 

is an exact solution with period 2p + 2 (p > 2) provided Eq. (|137|) is satisfied. From here it follows that such a 
solution is valid provided A 2 — —2A 3 — —2A 4 . 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 2p + 2) by 
randomly adding "a, a" between "0" and "a" or "a" and "0" , is an exact solution to Eq. {1} satisfying the constraint 
© provided Eq. (fnTjl is satisfied. 

(xxxii) tenth solution with period 6: <p = (..., a, —a, a, a, 0, a, ...) 
This is an exact solution to Eq. ([T]) satisfying the constraint ([2]) provided 

A - 4 = h 2 a 2 {A 3 + A 4 + A x - A 2 - A 5 - A 6 ) , (A - 2) = h 2 a 2 {A l +A 3 ~ A 4 ) , 

2(A-l) = h 2 a 2 {2A 1 +A 3 + A 2 + A e ), 2 = h 2 a 2 (Ar + A 6 ) . (138) 

From here it follows that such a solution is valid provided A 3 + A 5 = 2A 2 = 4A 4 . 
One can in fact generalize this solution and show that even 
(f> = (..., a, —a, a, (a, —a, a, p times), a, 0, a, ...) 

is an exact solution with period 3p + 3 (p > 1) provided Eq. (|138[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period Zp + 3) by 
randomly adding "—a" between two "a" , is an exact solution to Eq. (fTJ) satisfying the constraint © provided Eq. 
(|138p is satisfied. 

(xxxiii) eleventh solution with period 6: <fi = (..., a, 0, a, 0, —a, —a, ...) 
This is an exact solution to Eq. ([T|) satisfying the constraint © provided 

2(A - 3) = h 2 a 2 (2A 1 + A 3 + A 2 + A 6 ) , (A - 2) = h 2 a 2 (A 4 + A 3 - A 4 ) , 

2(A-1) = h 2 a 2 (2A 1 + A 3 + A 2 + A 6 ) , 2 = h 2 a 2 (A 5 + A 6 ) , A - 2 = h 2 a 2 A 1 . (139) 

From here it follows that such a solution is valid provided A 2 = A$ , = A 3 = A 4 = 0. 
One can in fact generalize this solution and show that even 
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(f> = (..., a, 0, (a, 0, p times), —a,— a,...) 

is an exact solution with period 2p + 2 (p > 2) provided Eq. (|139[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 2p + 2) by 
randomly adding "—a" between "0" and "a" or "a" and "0", is an exact solution to Eq. |T]) satisfying the constraint 
© provided Eq. (fl39|) is satisfied. 

(xxxiv) Solution with period 7: </> = (..., a, a, 0, —a, —a, 0, a, ...) 

This is an exact solution to Eq. (JTJ) satisfying the constraint ^ provided 

2(A- 1) = h 2 (2A 1 + A 3 + A 2 + A 6 ), a 2 = 1 . (140) 

One can in fact generalize this solution and show that even 
(f> = (..., a, a, 0, — a, — a, 0, (a, a, 0, —a, —a, 0, p times), a,...) 
is an exact solution with period 6p+ 1 (p > 1) provided Eq. (|140[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 6p + 1) by 
randomly adding as many number of "a" ("—a") as one wants between two "a" ("—a"), is an exact solution to Eq. 
(P} satisfying the constraint ((2) provided Eq. (|140[) is satisfied. 

(xxxv) second solution with period 7: <fi — (..., a, 0, — a, a, 0, — a, 0, ...) 
This is an exact solution to Eq. ([1]) satisfying the constraint ([2]) provided 

2(A-3) = h 2 a 2 {2A 1 + A 3 ~ A 2 ~ A 6 ) , A - 2 = h 2 a 2 A 1 . (141) 

One can in fact generalize this solution and show that even 
<f> = (..., a, 0, —a, a, o, —a, (a,0, — a, p times), 0,...) 

is an exact solution with period 3p + 1 (p > 2) provided Eq. (|141[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 3p+ 1) by 
randomly adding "0" between "—a" and "a" or "a" and "—a", is an exact solution to Eq. |T]) satisfying the constraint 
© provided Eq. ([Hi]) is satisfied. 

(xxxvi) third solution with period 7: <fi = (■■■, a, a> 0, —a, —a, a, —a, ...) 
This is an exact solution to Eq. (Q]) satisfying the constraint ((2j) provided 

2(A-l) = h 2 a 2 (2A 1 + A 3 + A 2 + A 6 ), A - 2 = h 2 a 2 A 1 , A - 4 = /i 2 a 2 (A 3 + A 4 + A 1 - A 2 - A 5 - A e ) . (142) 

From here it follows that such a solution is valid provided 2 A3 + A4 = A§ . 
One can in fact generalize this solution and show that even 
4> = (..., a, a, 0, —a, —a (a, a, 0, —a, —a, p times), a, —a, ...) 
is an exact solution with period 5p+ 2 (p > 1) provided Eq. (|142p is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 5p + 2) 
by randomly adding "—a, a" between "0" and "—a" or "a" and "0", is an exact solution to Eq. ([T]) satisfying the 
constraint © provided Eq. (|142p is satisfied. 

(xxxvii) fourth solution with period 7: <fi = (..., a, —a, 0, a, a, 0, —a, ...) 
This is an exact solution to Eq. (Q]) satisfying the constraint © provided 

2(A-1) = h 2 a 2 (2A 1 + A 3 +A 2 + A e ) , 2(A-3) = h 2 a 2 (2A 1 + A 3 -A 2 ~A 6 ) , A-4 - ?i 2 a 2 (^3+^4+^i-^2-A 5 -^ 6 ) . 

(143) 

From here it follows that such a solution is valid provided A 3 + 2A4 — 2A 5 . 
One can in fact generalize this solution and show that even 
(f> = (..., a, — a, 0, (a, —a, 0, p times), a, a, 0, — a, ...) 

is an exact solution with period 3p + 4 (p > 1) provided Eq. (|143[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 3p + 4) 
by randomly adding "a, —a" between "0" and "a" or "—a" and "0", is an exact solution to Eq. (JTJ) satisfying the 
constraint ([2]) provided Eq. (|143p is satisfied. 

(xxxviii) fifth solution with period 7: <fi = (..., a, a, 0, a, a, 0, —a, ...) 
This is an exact solution to Eq. (Q]) satisfying the constraint © provided 

2(A - 1) - h 2 a 2 (2A 1 +A 3 +A 2 + A 6 ) , 2(A - 3) = h 2 a 2 (2A 1 + A 3 - A 2 - A 6 ) , 

A- 2 = hV(A 3 + ,4i -A±), 2 = h 2 a 2 (A 5 +A 6 ). (144) 

From here it follows that such a solution is valid provided A 3 = 2A 4 , A 2 = A 5 . 
One can in fact generalize this solution and show that even 
4> =(..., a, a, 0, (a, a, 0, p times), —a,...) 
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is an exact solution with period ip + 1 (p > 2) provided Eq. (|144[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 3p + 1) by 
randomly adding "—a" between "0" and "a" or "a" and "0", is an exact solution to Eq. |T]) satisfying the constraint 
© provided Eq. (| 144ft is satisfied. 

(xxxix) sixth solution with period 7: <fi = (..., a, a, —a, a, a, —a, 0, ...) 
This is an exact solution to Eq. (p} satisfying the constraint ((2]) provided 

2(A - 1) = h 2 a 2 (2A 1 +A 3 + A 2 + A 6 ) , 2(A - 3) = h 2 a 2 (2A 1 + A 3 - A 2 - A 6 ) , 

A-2 = h 2 a 2 (A 3 + A 1 -A i ), A - 4 = h 2 a 2 (A 3 + A A + A 1 - A 2 - A 5 - A 6 ) . (145) 

From here it follows that such a solution is valid provided A 3 = 2A4 = A 5 . 
One can in fact generalize this solution and show that even 
4> =(..., a, a, —a, (a, a, —a, p times), 0,...) 

is an exact solution with period 3p + 1 (p > 2) provided Eq. (I145|) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 3p+ 1) by 
randomly adding "0" between "—a" and "a" or "a" and "—a" such that no two "0" are ever nearest or next-to-nearest 
neighbors, is an exact solution to Eq. (JTJ) satisfying the constraint ([2]) provided Eq. (| 145[) is satisfied. 

(xxxx) seventh solution with period 7: <fi = (..., a, 0, —a, a, a, a, —a, ...) 
This is an exact solution to Eq. ([T]) satisfying the constraint provided 

a 2 = l, 2(A-3) = h 2 (2A 1 +A 3 -A 2 -A 6 ), 

A-2 = h 2 {A 3 + A 1 -A 4 ), A-4 = h 2 (A 3 + A 4 + A 1 -A 2 -A 5 -A 6 ). (146) 

In view of the constraint {2]), such a solution is valid only if A 4 = 0, A 3 — A5. 
One can in fact generalize this solution and show that even 
<j> = (..., a, 0, —a, (a, 0, —a, p times), a, a, a, — a, ...) 

is an exact solution with period 3p + 4 (p > 1) provided Eq. (|146[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 3p + 4) 
by randomly adding "—a, a" between "a" and "0" or "0" and "—a", is an exact solution to Eq. {T]) satisfying the 
constraint j2|) provided Eq. (|146|) is satisfied. 

(xxxxi) eighth solution with period 7: <fi = (..., a, a, 0, a, 0, —a, —a, ...) 
This is an exact solution to Eq. (JT|) satisfying the constraint ([2]) provided 

A - 2 = h 2 a 2 A 1 , 2(A - 1) = h 2 a 2 {2A 1 +A 3 + A 2 + A 6 ) , 

A — 2 = h 2 a 2 (A 3 + Ai — A4) , 2 = /iV(A 5 + A 6 ) . (147) 

From here it follows that such a solution is valid provided A 3 = A4 , A$ = A 2 + A 3 . 
One can in fact generalize this solution and show that even 
4> =(..., a, a, 0, a, 0, (a, a, 0, a, 0, p times), —a, —a,...) 

is an exact solution with period 5p+ 2 (p > 1) provided Eq. (|147[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 5p + 2) by 
randomly adding "0" between any two "a" , is an exact solution to Eq. (fT]) satisfying the constraint provided Eq. 
(fT47j> is satisfied. 

(xxxxii) ninth solution with period 7: </> = (..., a, a, — a, 0, — a, 0, — a, ...) 
This is an exact solution to Eq. satisfying the constraint |2]) provided 

A - 2 = h 2 A x , 2(A - 3) = h 2 {2A 1 + A 3 - A 2 ~ A 6 ) , 

A-2 = h 2 (A 3 +A 1 -A 4 ), 2 = h 2 (A 5 +A 6 ). (148) 

In view of the constraint such a solution is valid only if A 2 = —2A 3 = —2A 4 , A5 = —3 A 3 . 
One can in fact generalize this solution and show that even 
4> =(..., a, a, —a, 0, —a, (a, a, —a, 0, —a, p times), 0,— a, ...) 
is an exact solution with period 5p+ 2 (p> 1) provided Eq. (1148|) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 5p + 2) by 
randomly adding "0" between any two "a" , is an exact solution to Eq. {1} satisfying the constraint ([2]) provided Eq. 
(|T48|) is satisfied. 

(xxxxiii) tenth solution with period 7: <\> = (..., a, —a, a, a, a, 0, a, ...) 
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This is an exact solution to Eq. (JT|) satisfying the constraint Q provided 

a 2 = l, A-4 = h 2 (A 1 +A 3 + A i -A 2 -A 5 -Ae), 2(A - 1) = h 2 (2A 1 + A 3 + A 2 + As) , 

A-2 = h 2 (A 3 +A 1 -A i ), 2 = h 2 (A 5 +A 6 ). (149) 

In view of the constraint such a solution is valid only if A 4 = A 2 = 0, A3 = — As- 
One can in fact generalize this solution and show that even 
<f> = (..., a, —a, a, a, a, 0, (a, —a, a, a, a, 0, p times), a,...) 

is an exact solution with period 6p+ 1 (p > 1) provided Eq. (|14£>[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 5p + 2) by 
randomly adding as many "a" as one wants between two "a" , is an exact solution to Eq. {TJ satisfying the constraint 
(121) provided Eq. (fT49|l is satisfied. 

(xxxxiv) eleventh solution with period 7: <\> = (..., a, 0, a, 0, a, a, —a, ...) 
This is an exact solution to Eq. (fT]) satisfying the constraint ([2]) provided 

A - 2 = h 2 a 2 A 1 , A - 4 = h 2 a 2 (A x + A 3 + A 4 - A 2 - A 5 - A3) , 2(A - 1) = h 2 a 2 {2A 1 + A 3 + A 2 + A 6 ) , 
A- 2 = h 2 a 2 (A 3 + A 1 - At) , 2(A- 3) = h 2 a 2 (2A 1 + A3 - A 2 - A3) , 2 = h 2 a 2 (A 5 + A3) • (150) 

From here it follows that such a solution is valid provided A2 = A3 = At = A5 = 0. 
One can in fact generalize this solution and show that even 
4> =(..., a, 0, (a, 0, p times), a, a, —a,...) 

is an exact solution with period 2p + 3 (p > 2) provided Eq. p50[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 2p + 3) by 
randomly adding "—a" between "a" and "0" or "0" and "a", is an exact solution to Eq. (fT]) satisfying the constraint 
© provided Eq. (Tr5T)|) is satisfied. 

(xxxxv) Solution with period 8: <f> = (..., a, —a, 0, a, a, a, 0, —a, ...) 
This is an exact solution to Eq. (TT]) satisfying the constraint (J5J) provided 

A - 4 = /i 2 (Ai + A3 + At - A 2 ~ A 5 - A3) , 2(A - 1) = h 2 {2A 1 + A 3 + A 2 + A3) , 

a 2 = l, 2(A- 3) = /i 2 (2Ai + A3 - A2 - As)- (151) 

In view of the constraint |2]), such a solution is valid only if As = 0, A3 = —2A4. 
One can in fact generalize this solution and show that even 
4> =(..., a, —a, 0, (a, —a, 0, p times), a, a, a, 0, — a, ...) 

is an exact solution with period 3p+ 5 (p > 1) provided Eq. (|151[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 3p + 5) 
by randomly adding "a,— a" between "—a" and "0" or "0" and "a", is an exact solution to Eq. ([T]) satisfying the 
constraint ^ provided Eq. (|151[) is satisfied. 

(xxxxvi) second solution with period 8: <\> = (..., a, 0, a, 0, a, a, a, —a, ...) 
This is an exact solution to Eq. ([T]) satisfying the constraint ([2]) provided 

A-4 = / l 2 (A 1 +A3+A 1 -A 2 -A3- A3), A-2 = h 2 A 1 , 2(A - 1) = h 2 (2A 1 + A 3 + A 2 + A 6 ) , 

a 2 = l, 2 = ^5+ A3), A-2 = h 2 (A 1 +A 3 -A 4 ), 2(A - 3) = h 2 (2A 1 + A 3 - A 2 - A3) . (152) 

In view of the constraint ([2]), such a solution is valid only if A 2 = A 3 = A4 — As = 0. 
One can in fact generalize this solution and show that even 
4> = (..., a, 0, a, 0, a, a, a, (a, 0, a, 0, a, a, a, p times), —a,...) 
is an exact solution with period 7p+ 1 (p > 1) provided Eq. (|152[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 7p+ 1) by 
randomly adding any combination of "—a" , "0" and "a" at any place, is an exact solution to Eq. (fT]) satisfying the 
constraint @ provided Eq. (|152[) is satisfied. 

(xxxxvii) Solution with period 9:0= (..., a, a, 0, —a, —a, 0, a, —a, 0, ...) 
This is an exact solution to Eq. (fT]) satisfying the constraint provided 

2(A-l) = / l V(2A+A3 + A 2 +A3), 2(A-3) = / l V(2A+A3-A 2 -A 3 ). (153) 

One can in fact generalize this solution and show that even 

4> = (•••, a, a, 0, —a, —a, 0, (a, a, 0, —a, —a, 0, p times), a, — a, 0,...) 

is an exact solution with period 6p + 3 (p > 1) provided Eq. (|153[) is satisfied. 



23 



One can also show that an aperiodic solution constructed from the above periodic solution (with period 6p + 3) 
by randomly adding "a, — a, 0" between "0" and "—a", is an exact solution to Eq. (Q]) satisfying the constraint © 
provided Eq. (|153|) is satisfied. 

(xxxxviii) second solution with period 9: 4> = (..., a, —a, 0, a, 0, a, a, 0, —a, ...) 
This is an exact solution to Eq. ([I]) satisfying the constraint ([2]) provided 

A - 4 = h 2 a 2 {A 1 + A 3 + A A - A 2 - A 5 - A 6 ) , A - 2 = a 2 h 2 A 1 , 2(A - 1) = h 2 a 2 {2A 1 + A 3 + A 2 + A 6 ) , 
2 = h 2 a 2 (A 5 +A 6 ), 2(A - 3) = h 2 a 2 (2 A 1 + A 3 - A 2 - A 6 ) . (154) 

From here it follows that such a solution is valid provided A 2 = A 4 = A 5 , A 3 = 0. 
One can in fact generalize this solution and show that even 
4> = (..., a, — a, 0, a, 0, a, a, 0, (a, — a, 0, a, 0, a, a, 0, p times), —a,...) 
is an exact solution with period 8p + 1 (p > 1) provided Eq. (|154[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 8p + 1) by 
randomly adding "0" between "a" and "—a", is an exact solution to Eq. fl} satisfying the constraint ([2]) provided 
Eq. (|154p is satisfied. 

(xxxxix) Solution with period 10: <f> = (..., a, a, 0, —a, —a. 0, a. —a, 0, a, ...) 
This is an exact solution to Eq. ([1]) satisfying the constraint ([2]) provided 

o a = l, 2(A-l) = / l 2 (2A 1 + A 3 + A 2 +A 6 ), 2(A-3) = /i 2 (2A 1 +A 3 -A 2 -A 6 ). (155) 

In view of the constraint |(3J|, such a solution is valid only if A 3 + 2A4 + 2A$ = 0. 
One can in fact generalize this solution and show that even 
4> = (..., a, a, 0, —a, —a, 0, a, —a, 0, (a, a, 0, —a, —a, 0, a, —a, 0, p times), a, ...) 
is an exact solution with period 9p+ 1 (p > 1) provided Eq. (|155[) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 9p+l) by 
randomly adding as many "a" as one wants between two "a" and by randomly adding as many "—a" between two 
"—a" is an exact solution to Eq. fTJ) satisfying the constraint @ provided Eq. (|155|) is satisfied. 

(xxxxx) second solution with period 10: <fr = (..., a, —a, 0, a, 0, a, a, a, 0, —a, ...) 

This is an exact solution to Eq. (|TJ) satisfying the constraint {2]| provided 

A-4 = /i 2 (^ 1 +A 3 +A4-A 2 -A 5 -A 6 ), A-2 = /i 2 A 1; 2(A - 1) = h 2 (2A 1 + A 3 + A 2 + A 6 ) , 

2 = h 2 {A 5 +A e ), a 2 = l, 2(A-3) = /i 2 (2 J 4 1 +A 3 -A 2 -yl 6 ). (156) 

In view of the constraint ((2]), such a solution is valid only if A2 = A3 = A4 = A5 = 0. 
One can in fact generalize this solution and show that even 
4> =(..., a, —a, 0, a, 0, a, a, a, 0, (a, — a, 0, a, 0, a, a, a, 0, p times), —a,...) 
is an exact solution with period 9p+ 1 (p > 1) provided Eq. (| 1 56|) is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 9p+ 1) by 
randomly adding as many "a" as one wants between two "a" and by randomly adding "0" between "—a" and "a" or 
"a" and "—a" is an exact solution to Eq. |T]) satisfying the constraint ^ provided Eq. (|156[) is satisfied. 

(xxxxxi) Solution with period 11: <f) = (■■-, a , — a , 0, a, a, a, 0, —a, a, a, —a, ...) 

This is an exact solution to Eq. (JTJ) satisfying the constraint ^ provided 

A - 4 = h 2 {A x + A 3 + A 4 - A 2 ~ A 5 - A 6 ) , 2(A - 1) = h 2 (2A 1 + A 3 + A 2 + A 6 ) , 

a 2 = l, 2(A- 3) = h 2 (2A 1 + A 3 - A 2 - A 6 ) , A - 2 = h 2 {A 1 + A 3 - A 4 ) . (157) 

In view of the constraint such a solution is valid only if A 3 = A4 = A$ = 0. 
One can in fact generalize this solution and show that even 
<j> = (..., a, —a, 0, a, a, a, 0, —a, a, a, (a, —a, 0, a, a, a, 0, —a, a, a, p times), — a, ...) 
is an exact solution with period 10p+ 1 (p > 1) provided Eq. (|157p is satisfied. 

One can also show that an aperiodic solution constructed from the above periodic solution (with period 10p+ 1) by 
randomly adding as many "a" as one wants between two "a" , is an exact solution to Eq. {1} satisfying the constraint 
© provided Eq. (fH)T)) is satisfied. 
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V. TWO-POINT MAPS 



A. General maps that include integration constant 

As pointed out in Sec. QJ while Model 2, for any value of 5 and 7, has the first integral with the integration constant 
C expressed by Eq. for the Speight and Ward Model 7, the known two point map Eq. (fTTj) does not contain the 
integration constant C . 

Remarkably, in case only A 2 and A4 are nonzero (combination of Models 3 and 6, i.e. an admixture of Bender- Tovbis 
and BOP nonlinearities), the discrete model Eq. ([1]) has the following first integral with the integration constant C 



where 



K + x = (2 - A) ^:^^^ 1 , f(K) = ~ Xfi + 4 J, (158) 

z (2-A) 2 -C/A4t _ CY 2 + (2 - A) 2 (l - Z 2 ) Y _ h * (A , AZ) (im 
Z - 2(2- A) + Ch^A 2 A^ X - (2~~A)Y ' Y-h{A A + A 3 Z). (159) 

As expected, in the limit of either A 2 — or A4 — 0, the nonlinear map reduces to the maps derived earlier in [l2l.[l4l]. 
Thus the nonlinear map defined by Eqs. (|158[) and (|159[) generalizes the maps considered in [12, [H| and it allows one 
to construct static solutions to Eq. ([I]) iteratively in a way similar to those studies. 
Remarkably, in case A ^ 2, this first integral can be further factorized as 

W(<t> n , <£„+i)(2 - A) (2 - A - Y4>1) = , (160) 

where 

Y CY 

W(<p n , <j> n+1 ) = + <g +1 - — ^4>14>1 +1 - - ■ (161) 

Since the vanishing of the third bracket in Eq. (|160p is a trivial possibility, for A 7^ 2, effectively the first integral in 
this case acquires rather simple form, W (4> n , 0n+i) = 0, which is precisely of the Quispel form as given in [22j . As 
expected, in the special case when A4 = and hence h 2 A 2 = A, the first integral W((j) n , <j> n +i) — agrees with that 
obtained from the map Eq. ^ in case 7 = S = 0. We also note that Eq. ([1]) with only A 2 and A4 nonzero can be 
expressed in terms of Eq. (|161[) as follows 

= nv,A~ k A T {W^n, <f>n+l) ~ W{<t> n -i, <j> n ) + [(& +1 W(</»„-l, <f>n) ~ <f>n+l)] } ■ (162) 

2Z(0 n+ i - (pn-i) L 2-A L J J 

It is now clear that the static solutions to Eq. ([1]) can be found from the two-point map W(<p n , 4> n +i) = 0. Note 
that for small h one has Z « (2 — A)/2 and the last term in the curly bracket of Eq. (|162[) can be neglected and one 
obtains the equation similar to Eq. of Model 2. In other words, the model with only A 2 and A4 nonzero can be 
regarded as the Model 2 modified by the O-term (the last term in the curly bracket), i.e. the term which disappears 
in the continuum limit and vanishes upon substituting W(<fi n , 4> n +i) = [12j. 

We have checked that from the two-point map W(4> n , 4> n +i) — one can obtain the staggered as well as the 
nonstaggered JEF solutions dn, cn and sn derived in Sec. Ill from the three-point static Eq. ([T]) in case only A 2 and 
A4 are nonzero. The following local identities are helpful in these derivations [23| 

dn 2 (x, m)dn 2 (x+a, m) + cs 2 (a)[dn 2 (a;, m) + dn 2 (a; + a, m)] — 2ds(a, m)ns(a, m)dn(ir, m)dn(a; + a, m) + l — m = , (163) 

m 2 cn 2 (x, m)cn 2 (x+a, m)+mds 2 (a)[cn 2 (x, m)+cn 2 (x+a, m)] — 2mcs(a, m)ns(a, m)cn(x, m)cn(x+a, to)— m(l— m) = , 

(164) 

TO 2 sn 2 (x, TO.)sn 2 (x+a, m)— TOns 2 (a)[sn 2 (x, TO,)+sn 2 (x+a, m)]+2TOCs(a, m)ds(a, m)sn(x, m)sn(x+a, m)+m = . (165) 

It is worth emphasizing here that as demonstrated in [l4| , first integral can always be constructed from the known 
JEF solutions. In the continuum limit, the DFI Eq. (|158[) and the DFI W(4>n, 4>n+i) — as given by Eq. (|161[) . 
reduce to Eq. which is the first integral of the static continuum <fi 4 equation. Relation between the solutions 
obtained from the map W(4> n ,4>n+i) = iteratively and JEF solutions reported in Sec. IIIIBI can be established as 
follows. JEF solutions include two integration constants, the arbitrary shift xq and the modulus to. In the map 
W(4> n , 4> n +i) = 0, the role of to is played by C and the role of xo is played by the initial value of the map <pQ. Relation 
between m and C can be found if one obtains the JEF solutions not only from the three-point static problem of Eq. 
([1]) but also from the two-point problem W(4> n ,4>n+i) = fl2l. [l4T|. 
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B. Two-point maps for the cases (iv) to (vii) 

What about a universal map for the cases (iv) to (vii)? Unfortunately, so far we have not been able to find one for 
any of these cases. However, as we show now, corresponding to any static JEF solution, we can always generate a 
map. In particular, we now obtain maps corresponding to the three JEF solutions sn, cn and dn and show that they 
can always be factorized and reduced to Quispel [22[ form. 



1. sn solution 



As shown in the last section, one of the exact JEF solution in cases (iv) to (vii) is 

4> n = AsD.[hf3(n + xq), m] , (166) 

with arbitrary xo and < to < 1. On using the well known identity 

sn(u, m) cn(v, m) dsx(v, m) + sn(v, m) cn(u, m) dn(it, m) 

sn(u + v,m) = ' y > > v L_!_^ (167) 

1 — m sn z (it, m) sn z [y, to) 

it immediately follows that 4> n +i and 4> n are related by the map 

(f> n cn(h/3, m)dn(h/3, to) + sn(hf3, m) y/A 2 - (1 + m)4% + m(<f>*/A 2 ) 



Vn+l 



(168) 



1 — msn 2 (/i/3, mjcfx^/A 2 
This can be simplified and put in the following factorized form 

[A 2 - msn 2 (hp, m)<f>l]W sn {<j> n , cj> n+l ) = , (169) 

where 

W an (4> n ,<f> n+ i) = A 2 (4> 2 n + <j) 2 n+1 ) - msn 2 (hp, m)4> 2 n (}? n+1 - 2A 2 cn(/i/3, m)dn(ft/3, m)(f> n <j) n + 1 - A 4 sn 2 (hf3 7 m) . (170) 

Since the vanishing of the first bracket in Eq. (|169[) is a trivial possibility, effectively the map in this case is given by 
W r sra (0n, <t>n+i) = which is precisely of the Quispel form as given in [22| . 

For any of the cases (iv) to (vii) , the sn solution Eq. (|166p can be derived from this map for any initial value taken 
from —A < (j)Q < A. In this map, < m < 1 is the integration constant. For given h, parameters of the map j3 and 
A are related to the model parameters. 



2. dn solution 



Similarly, on using the identity 

, , . dn(u, m) dn(u, m) — msn(v, m) + cn(v, m) sn(u, m) cn(tt, m) . _ . 

dn(w + V,?7l) = J? N 2? ^ ' ( 171 ) 

1 — to sn z (u, to) sn z (w, to) 

we find the map for the dn solution 

(p n = Adii[h(3(n + x ),m] , (172) 

in the form 

0ndn(fy3, to) - sn(h/3, m)cn(/i/3, m)^A 2 {m - 1) + (2 - m)<j) 2 n - {$JA?) 



(173) 



1 - sn 2 (/i/3,TO)(l -<j>l/A 2 ) 
This can be simplified and put in the following factorized form 

[A 2 dn 2 (h(3,m) + tosu 2 ^, to) ( /. 2 ]^„(0„, 0„+i) =0, (174) 

where 

W d n(4>n, 4>n+i) = A 2 dn 2 (h(3, to) {<j) 2 n + 4>n + i) + msn 2 (hj3, m)<)) 2 n (f> 2 n+1 - 2A 2 cn(hf3,m)cf) n (f) n+1 - A 4 (l - m)sn 2 (hf3,m) . 

. (175) 

Since the vanishing of the first bracket in Eq. (|174|) is a trivial possibility, effectively the map in this case is given by 
Wdni&n, 0n+i ) = which is precisely of the Quispel form as given in [22| |. 
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3. cn solution 

Similarly, on using the identity 

, . enfu, m) cn(v, m) — sn(v, m) + dn(v, m) sa(u, m) dn(n, m) ,«„„ s 

cn(u + v,m) = 1 ; v > y ' > \\ ' v ; v ; , (176) 

1 — m sn z (u, m) sn z (-u, m) 

we find the map for the cn solution 

ct> n = Acn[hf3(n + x Q ),m] , (177) 

in the form 



<M+1 



<f> n cn(h0,m) - sn{hf3,m)dn(h/3,m) v / A 2 (l - m) + (2m - 1)02 _ m (cj) 4 n /A 2 ) 



(178) 



1 - msn 2 (hf3,m)(l - 4> 2 /A 2 ) 
This can be simplified and put in the following factorized form 

[A 2 cn 2 (hf3,m) + sn 2 {h(3, m)0 2 ]VF c „(0„, 0„+i) = 0, (179) 

where 

W cn {<j> n ,<j> n+1 ) = A 2 cn 2 {h(3,m){<l> 2 n +(g + {)+sn 2 {h(3,m)^ . (180) 

Since the vanishing of the first bracket in Eq. (|179p is a trivial possibility, effectively the map in this case is given by 
W cn (<f) n , 4>n+i) — which is precisely of the Quispel form as given in [221 ]. 

Summarizing, unlike the cases (i) to (iii), we do not have a universal map for the cases (iv) to (vii). However, in 
these four cases, one can always obtain a map by starting from sn, cn or dn solutions and remarkably, the maps in all 
three cases is effectively of the Quispel form. Obviously, one can also obtain the staggered as well as the nonstaggered 
JEF solution from the relevant two-point map. 

It is clear from these arguments that for any discrete model that admits a JEF solution, one can easily construct 
the corresponding two-point map from which that JEF solution can be iteratively generated [l4j |. Further, such a map 
is effectively of the Quispel form. However, the problem of finding a universal map from which any static solution can 
be generated is not trivial even when the JEF solutions are known (note that not all TI models admit JEF solutions) . 

C. Particular factorized static problems 

In some cases, the static problem can be factorized and one can obtain some of the exact solutions, such as those 
presented in Sec. IIVI from this lower order algebraic equation. As an illustration, we now present a few examples of 
such factorized static problems. 

As an example, we note that the two-point map for the Model 2 with 7 = and 5 = 1/4 is given by Eq. ([6]) 



{<t>n+l ~ <t>nf + A0„0n+1 - ^0„<?Wl (4>n + 4>n+l) ~ \ + C = ■ ( :l ,V| ! 

If the integration constant is chosen as Ch 2 = (A — 4) 2 /(2A), then Eq. (|181[) factorizes as 



4 ( A , 

A V 1 ~ l^ 71 ^ 1 



A - 2 - j (<j> 2 n + 4>l +1 ) 



0. (182) 



The last multiplier of Eq. (|182p generates the TI four-periodic solution <j> n = (..., a, b, —a, —6, ...) with A(a 2 + b 2 ) = 
4(A — 2). This solution can also be written in the form of the four-periodic TI sine solution as given by Eq. (f8"2"| and 
Eq. Further, the last multiplier also generates the staggered TI sine solution with h(3 — n/2, although it can be 

noted that for the four-periodic solution, non-staggered and staggered forms coincide after n is substituted with —n. 

As a second example we take the Bender- Tovbis model, i.e. the Model 2 with 7 = 5 = 0. In this case, the two-point 
map Eq. © with A = 2 and C = 1 is factorized as 

(l-0 2 l _ 1 )(0'-l) = O. (183) 



This example explains how one can obtain aperiodic solutions from factorized maps. Indeed, any sequence of ±1 
satisfies Eq. (|183[) and hence the Bender- Tovbis model at A = 2 and C = 1. To obtain this aperiodic solution one 
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can use either of the multipliers of Eq. (|183[) . This is possible because the solutions are derived from the two-point 
rather than the three-point map. 

In the above examples the short-periodic solutions were obtained from the factorized two-point map. More examples 
of this sort can be found in [l4| where short-periodic solutions are derived from the two-point map Eq. (|158[) factorized 
for A = 2 for the cases A 2 = A and A4 = A. 

As another example, the Speight and Ward Model 7 can be written in the form [12] 

d d 

4>n = -«(0 n -l, <l>n)-^l-v{4>n-l, <f>n) ~ K0n; 0n+l)^— u (0n, 0n+l) , (184) 

where 

/, , x 0„-0„_l 1 . (4>l-l + 0n-10n + 4>l) „2 6A / 1QK s 

«(^-i.^)=— ff — 71+ ^ > ff= 6^A' (185) 

from which it is clear that the equation u(0 n _i, <j) n ) — generates static kink and inverted kink solutions of their 
model. 

The following two examples are interesting because they give the short-periodic solutions, not from a two-point 
map, but, from a set of two finite-difference equations. 

We note that the Speight and Ward Model 7 can be written in the form 



071 = 

~ Jg(^-1 + %M(&-1 + 071-107, + <t>l) ~ + Wn){4>l + 0n0n+l + f n +l)- ( 186 ) 

If the following two-point equation holds 

+ K-iK + 4>l = 6 ^x^' ( 187 ) 

then the static version of Eq. (|186p reduces to 

(A - 6) (0 n _i + <f>n+ <l>n+i) = . (188) 

Two-point map Eq. (| 1 8T[) generates the exact three-periodic sin as well as the six-periodic staggered sin solutions to 
the Speight and Ward model [see Eqs. (|8"2"]). (f9"T)) and (|9"8)) ] and for these solutions Eq. (|188[) is also satisfied. Further, 
Eqs. (|187| and l|188|) for the case of A = 6 also generate the short periodic solutions n = (...,a,0, ..), (...,a, — a,0, ...), 
(..., a, —a, a, — a, 0, ...), (..., a, 0, a, 0, —a, ...), (..., a, —a, 0, — a, a, 0, ...), (a, 0, a, 0, a, —a), and (a, 0, —a, a, 0, —a, 0). 
Similarly, the equation of motion for the Hamiltonian Model 10 can be written in the form 

0n = Tp (<t>n+l + 0n-l - 20n ) + A0„ 

-A(aa0 n -i + 2ai0„)(0 r ^_ 1 + — 0„_ 1 0„ + <p\) - A(a 2 0,!+i + 2ai0 n )(0L_ 1 + — 0„ +1 0„ + 0^) . (189) 

ai ai 



In case the following two-point equation holds 



d + — 0r!-10r! +4>l=H, (190) 
Oil 



then the static version of Eq. (|189p reduces to 



provided 



n+ i)(l-Aifa 2 ) + (A-2-4aiAif)0„=Q, H = ±- 2 a2 J. a \ , (191) 

A(4ai — Q^joLx) 



a 3 = at + —2- . (192) 
2ai 



In the special case of a,\ — a 2 — 1/18, this reduces to the SW model. The two-point map Eq. (| 1 90[) generates the 
exact sine as well as the staggered sine solutions to the Hamiltonian Model 10 with cos(h(3) = ^fa 2 /(2ai), and for 
these solutions Eq. (|19ip is also satisfied. 
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Once again, it is interesting that for the factorization Eq. (|187[) and Eq. (|188[) . and also for the factorization Eq. 
(|190[) and Eq. (|191[) . one has to satisfy two lower-order finite-difference equations simultaneously, and one of those 
equations is a two-point one while another is a three-point one. 

All factorized problems discussed in this section do not contain the integration constant and thus they generate 
only particular solutions. Some of them are TI solutions, for example, the three-periodic solution derivable from Eq. 
(| 187)1 and Eq. (|188|) . while others are not, for example, arbitrary sequence of ±1, derivable from Eq. (|183|) . 

Solutions constructed in this section from factorized problems do not survive in the continuum limit because the 
factorized problems do not reduce to Eq. ([T]) or Eq. ([J} in the continuum limit. 



VI. GOLDSTONE MODES 



A. Goldstone mode of a TI static solution 



Let 4% be a static solution to Eq. ([I]). To study the dynamics in the vicinity of this solution we substitute the 
ansatz 4> n (t) — <^>„ + £ n (t) into Eq. ([T]), and obtain the following linearized equation 



(193) 
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Looking for solutions of Eq. ()193|) of the form e n (t) — U n exp(±iujt) we come to the eigen- value problem 

[K] U = -iu 2 U , 



^<_ 1 «_ 1 +2^ +1 )-^ £ (0° +1 ) 2 . 



(194) 



(195) 



where vector U contains U n and the nonzero coefficients of matrix [K] are given by Eq. (|194p . 

If is a TI static solution then it can be shifted along the chain by an arbitrary x , 0° = <fi(n + x ). The 
eigenvector corresponding to the zero- frequency translational Goldstone mode, Uq, has components 4>' n , where prime 
means derivative of </> with respect to its argument. To confirm that, we substitute Uq = {^'n} into Eq. (|195p with 



■jj- 



and obtain 



-1 H~ -^n,n^ 



(196) 



The last expression is an identity because it coincides with the derivative of static version of Eq. ([T]) with respect to 
xq, and such a derivation is possible for the TI static solution, which is an equilibrium solution for any xq. We thus 
have proved that any TI static solution has the zero- frequency translational mode Uq = {</4}- Particularly, for any 
static JEF, hyperbolic or trigonometric function solutions given in Sec. IIII1 one can easily find the corresponding TI 
mode as it is proportional to the derivative of the solution with respect to its argument. 

Looking for solutions of Eq. (|193[) and Eq. (|194[) with 0° = 1 of the form of small-amplitude phonons, e n {t) ~ 
exp(i/cn ± itot), where k denotes wavenumber and w is frequency, one obtains the spectrum of the vacuum for the 
discrete model of Eq. (JTJ) , 



2A + 2 



h 2 



- A 2 - 2A 3 - 2Aa - 3A 5 - 3A 6 



sin 



(197) 



B. Goldstone modes of some short-period static solutions 



If a static solution does not possess the zero-frequency translational Goldstone mode, then this solution is not a 
TI one. The opposite, in general is not true, i.e., a particular static solution may have the Goldstone mode only at 
certain positions with respect to the lattice xq, but a TI solution must have such a mode at any xq. It is interesting 
to check whether the short-periodic solutions derived in Sec. [TV] can have the Goldstone mode. 



29 



Four periodic solution of the form a, b, —a, —b, as found in Sec. IIV[ exists in case A\ ^ under the constraint 
Ai — A 3 + A4 — and has a 2 = (A — 2)/{A 1 ti 2 ) — b 2 . This one-parameter solution is a TI solution because it can be 
expressed in the form </>„ = Acos[n(n + xq)/2] with arbitrary shift xo and A 2 = (A — 2) / (Aih 2 ) . Being a TI solution, 
it possesses the Goldstone mode at any xo, as it was demonstrated in Sec. IVI A[ Note that Model 2, Model 8, Model 
10 and Speight and Ward Model 7 have A\ ^ but the constraint A\ — A3 + A4 = is not satisfied for the SW model 
7, while it is satisfied for Model 2 (at arbitrary 7 and arbitrary nonzero 5), Model 8 (at 7 = 0, a = 1/4) and Model 
10 (at ax — 0:3). Thus, while Models 2,8 and 10 support this TI four-periodic solution, Model 7 does not. 

It may be noted here that in case A = 2 and A\ = Q,A 3 = A4, then the TI solution is a two-parameter solution 
which can be expressed in the form (f>„ = Acos[i:(n + x$)/2] with an arbitrary shift x$ and an arbitrary amplitude 
A. This is satisfied by Model 9 at fj = as well as by models where only A 3 , A4, ^5 or where only A 2 , A 3 , A4, A 5 are 
nonzero. 

Four periodic solution of the form a, 0, — a, 0, as found in Sec. HVi exists for A\ 7^ and has a 2 = 
(A — 2) / (Aih 2 ). Inserting this solution into Eq. (|195[) one finds 
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a = - 

7=^-^3 + ^)-^, 

Characteristic equation of Eq. (| 198(1 is 

(a + w 2 ) (7 + uj 2 ) [(a + uj 2 ) (7 + oj 2 ) - 4/5(5] = 0, 
and zero-frequency modes are possible when any of the following three conditions is satisfied 

a = 0, 7 = 0, cry - 4/3(5 = 0. 



(199) 



(200) 



(201) 



The considered four-periodic solution can also be expressed as </>„ = Acos(7rn/2) with A 2 = (A — 2)/(Aih 2 ). The 
expected Goldstone mode is proportional to the derivative of this function with respect to its argument, and thus, 
should have the form Uq = 11% = 0, Ui = — U3 = k with arbitrary k ^ 0. Substituting this into Eq. (|198|) with 
uj 2 = one finds that the Goldstone mode corresponds to 7 = 0. Taking into account Eq. (|199|) . we conclude that 
the considered four-periodic solution possesses the Goldstone mode under the condition 

A x 



A, + Aa 



0. 



(202) 



Particularly, the Speight and Ward Model 7 supports the considered four-periodic solution but it does not satisfy 
the condition Eq. (|202[) so that the solution is not a TI one. On the other hand, Model 2 supports the considered 
four-periodic solution and it satisfies the condition of Eq. (I202[) so that the solution possesses the Goldstone mode. 

It may be noted here that the solution (...,a,0,-a,0,...) also exists in case A = 2 and A\ = 0, and in that case a is 
any real number. In this case the solution can be expressed in the form </>„ = acos[7m/2]. Further, such a solution is 
always a TI one since according to Eq. (|201|) . 7 is proportional to A — 2 and hence is equal to 0. Note that such a 
solution exists in almost all the models at A = 2. For example it exists in Model 2 at 6 — 0, Model 3, Model 4, Model 
6, Model 8 (at a = or/and 7 = —(3 — 1/2), Model 9 and Model 10 (at a% = 0) as well as the five models discussed 
in Sec. III. 

Four periodic solution of the form ...,a,a,—a,—a,..., as found in Sec. HVi exists for A\ + A 3 — A4 7^ and 
has h 2 a 2 = (A - 2)/(Ai + A 3 - A4). The solution can also be presented as <j) n = v / 2^4cos(7rn/2 + tt/4) where 
A 2 = (A - 2)/[h 2 (A 1 + A 3 - A4)] , so that the Goldstone mode is Uo = —U\ = -U2 = U 3 = k with an arbitrary k ^ 0. 
Performing the calculations similar to the previous case one finds that the considered four-periodic solution possesses 
the Goldstone mode under the same condition of Eq. <|202|) . Again, this solution has the Goldstone mode in Model 2 
but not in the Speight and Ward Model 7, so that in the latter model the solution is not a TI one. 

It may be noted here that the solution (...,a,a,-a,-a,...) also exists in case A = 2 and A\ + A 3 — A4, and in this case 
a is any real number. In this case the solution can be expressed in the form <p n = \[2a cos[7rn/2 + 7r/4]. Further, such 
a solution is always a TI one since according to Eq. (|202|) , 7 is proportional to A — 2 and hence is equal to 0. It is 
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FIG. 1: Difference between (j>e and <f>\ as function of 0o for the static solutions constructed for the Speight and Ward Model 7 
from the three-point map (j) n +i = /(0n-i,0n) for chosen 4>q and numerically found <j)i such that 4>z, = (j>o- If the three-point 
problem is reducible to a two-point problem 0„+i = g(4> n ) then from having <j>5 — 4>o one must also have cj>$ = cj>i, which is not 
the case and thus, for this solution the two-point reduction is impossible. Static solutions generated by the three-point map 
for the points marked with a to e are shown in the corresponding panels of Fig. [2] 

easily checked that such a solution at A = 2 exists in Model 2 (at 6 = and arbitrary 7) , Model 8 (at 8a 2 = 2a + (3) , 
Model 9 (at = 0) and Model 10 (at a x = -a 3 ). 

Three periodic solution of the form a, 0, —a, as found in Sec. IIV1 exists for 2A\ — A 2 + A 3 — Aq ^ and has 
h 2 a 2 = 2(A — 3)/(2Ai — A2 + A 3 — Aq). Thus such a solution can be expressed by <p n = (2/>/3)A cos(27rn/3 — 7r/6) 
with A 2 = [2(A - 2)}/[h 2 {2A 1 - A 2 + A 3 - A 6 )}. The Goldstone mode is U = U 2 = k, U x = -p, and the ratio k/p 
can be found but it is not important for our analysis. Calculations similar to that for the four-periodic solutions give 
that the considered three-periodic solution possesses the Goldstone mode or, in other words, is a TI solution under 
the condition 

( - AA\ - A 2 + A 2 + 2A 2 - 2A 1 A 6 + 2A 1 A 2 + 4A X A 3 - iA x A 4 - A 2 A 6 + A 3 A 6 

+2A A A 6 - 3A 5 A 6 - A 2 A 3 + 2A 3 A 4 + A 2 A 5 - 2A|A 5 ) (A 2 - 3) 2 = . (203) 

It is easily checked that both Model 2 and Speight and Ward Model 7 support the three periodic solution and meet 
the condition Eq. (|203[) . Thus the three-periodic solution possesses the Goldstone mode in both Model 2 and Model 
7, but, as mentioned above, one cannot claim that this solution is a TI one in these models unless one demonstrates 
that the Goldstone mode exists for any shift of the three-periodic solution along the lattice. 

It may be noted here that the solution (...,a,0,-a,...) also exists in case A = 2 and 2A\ + A 3 = A 2 + Aq, and in this 
case a is any real number. In this case the solution can be expressed in the form <f> n = (2/\/3)acos[27rn/3 — 7r/6]. 
Further, such a solution is always a TI one since A = 3. Thus, such a solution exists at A = 2 in Model 2 (in case 
27 + 86 = 1), Model 6, Model 8 (in case Aa[/3 + 7] = 3a 2 + 7 2 + £7), Model 10 (in case 2a x + a 3 = 2a 2 ), and in 
Models where only A 2 , A 3 , A$ or only A 2 , A 3 , A4, A5 are nonzero. 

VII. NUMERICAL RESULTS 
A. Five-periodic static solution in Speight and Ward Model 7 

Here we address the problem of integrability of the Speight and Ward (SW) Model 7. As it was mentioned, this 
model has the two-point map Eq. (|llj) to derive the kink solution but a general two-point map, that includes the 
integration constant as a free parameter, for this model is not known. We will give numerical evidence that the Model 
7 is not integrable and the two-point map for obtaining a two-parameter set of static solutions cannot be constructed. 
For this purpose we will make an attempt to construct the static five-periodic solutions to this model and will observe 
that it can be constructed only for highly symmetric positions of the solution with respect to the lattice so that the 
solution is not a TI one and it possesses the Peierls-Nabarro potential. We have chosen the five-periodic solution for 
this study because it is relatively simple and it corresponds to a non-factorized static problem. Simple short-periodic 
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FIG. 2: Static solutions for the Speight and Ward Model 7 generated from the three-point map cj> n +i = f(cj> n -i,4>n) for chosen 
4>o and numerically found 4>i such that 0s = 4>q. Panels (a) to (e) show the results for the values of <f>o marked with a to e 
in Fig. [T] Solutions in (b) and (d) are not five-periodic solutions because 4>s differs from <j>i by the amount shown in Fig. [T] 
Highly symmetric solutions in (a), (c), and (e) are the five-periodic ones and for them (j>6 — <t>i, as it can be seen from Fig. [T] 
Solution in (a) [and in (e)] corresponds to a minimum of the Peierls-Nabarro potential while the one in (c) to a maximum of 
this potential. 

solutions described in Sec. IIVI are not suitable for this study because they are obtained from low-order algebraic 
equations, i.e., from factorized static problems that do not represent the model in its general formulation. 

We set for the model parameters A = 1, h = 1.3, and Ak corresponding to the Model 7 as given in the Introduction. 
The three-point static problem of Eq. (fTJ) is written in the form of the map 4> n +i — f(4>n—~i> 4>n) which, for given 
4>q and (f>\, generates a static solution. For chosen 4>q we numerically find 4>i such that in the itcratively obtained 
solution cf>5 — <fio and check whether in this case we also have 4>e — 4>i- If the three-point map (j) n +i = f{4>n-i 1 4'n) 
is reducible to a two-point map 4> n +\ — g{4> n ) then having <^5 = </>o we must also have <f>e = <f>\. However, as it can 
be seen from Fig. [TJ 6 — is equal to zero only for a discrete set of 4>q. In Fig. [2] (a) to (e) we plot the structures 
generated by the three-point map for different 4>q indicated in Fig. [T]by, correspondingly, letters a to e. Recall that 
for all structures presented in Fig. [5] we have ^5 = <fio but the condition tpg = <\>\ is fulfilled only for the highly 
symmetric structures shown in (a), (c), and (e) [structures in (a) and (e) are equivalent]. Solutions in (b) and (d) are 
not five-periodic solutions because 4>6 differs from (f>i by the amount shown in Fig. [TJ Solutions in (b) and (d) are 
modulated five-periodic structures but this cannot be seen in Fig. [2] because the period of the modulated structure is 
very large. 

Static solutions shown in Fig. [21(a) and (c) are the five-periodic equilibrium solutions for which the small-amplitude 
vibrational spectrum can be calculated as described in Sec. IVI Al Doing so we find that these five-periodic structures 
do not possess the zero-frequency Goldstone mode, but they have a nearly translational mode with frequency lu — 
0.0036 for the structure in Fig. [U(a) and purely imaginary frequency w — 0.0036i for the structure in Fig. d](c). Thus, 
the five-periodic structure in the Speight and Ward model 7 is not a TI one and it experiences the Peierls-Nabarro 
potential with a minimum energy corresponding to the structure in (a) and a maximum energy corresponding to the 
structure in (c). 

We have studied some other periodic solutions, for example, seven- and eight-periodic ones and have obtained the 
results qualitatively similar to that for the five-periodic structure. In all studied cases the structures had nonzero 
Peierls-Nabarro potential. 

We conclude that the static solutions supported by the Speight and Ward Model 7, except for the kinks, anti-kinks, 
sine and staggered-sine solutions, usually have the Peierls-Nabarro potential. The corresponding three-point static 
problem is non-integrable and cannot be reduced to a two-point problem, again, except for the kink solutions. 

This result is not surprising at all and it could be expected taking into account that the derivation of the two-point 
map Eq. (|11[) . from which the kink solution can be derived, was done for the integration constant C — in Eq. 
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FIG. 3: Static kink in the case (iii) with only A2 and A4, nonzero: (a) model parameters At and kink parameter (3 as functions 
of h; (b) on-site kink at h — 0.8; (c) inter-site kink at h = 0.8. This model admits TI static solutions, including the kink 
solution, at constant (h- independent) model parameters Ak, as can be seen in (a). 



The resulting discrete model supports the TI solutions only for this particular value of the integration constant, and 
those solutions are kinks. 

Hamiltonian TI discretization of the Klein-Gordon field that generalizes the Speight and Ward model has been 
derived in Sec. II C of 12]. The model includes the integration constant and thus supports a two-parameter set of 
TI static solutions, although it is rather complex even for the cubic nonlincarity. 

It is worth noting here that the Speight and Ward Model 7, apart from the five-periodic solutions shown in Fig. [2] 
(a), (c), and (e), supports the following non-TI five-periodic solutions: 



Solution (xiv) 


= 


(M,i,- 


1,-1) if A 


= 9/2; 


Solution (xix) 


<t> = 


(1,1,0,- 


1,0) if A = 


18/7; 


Solution (xvi) 


cj> = 


(a, —a, a, 


-a,0) if A 


= 6 and a 2 — 3 ; 


Solution (xvi) 


4> = 


(a, 0, a, 


- a) if A = 


6 and a 2 = 3 ; 


Solution (xvi) 


= 


(a, a, —a, 


-o,0) if A 


= 12 and a 2 = 3/2 



The solutions shown in Fig. [5] exist for a continuously varying A (at least within a range of A values) while the above 
five solutions exist for a fixed A. 



B. Static kinks 



Here, after a brief discussion on the JEF solutions, we focus on the analysis of the kink solutions because they are 
discussed in applications more often than the periodic solutions. 

As it was mentioned, the JEF solutions and their hyperbolic function limit solutions like kink and pulse exist in 
the model Eq. (J7XJ) m the seven cases, of which the first three cases with (i) only A2 nonzero, (ii) only A4 nonzero, and 
(iii) only Ai and A4 nonzero are qualitatively different from the other four cases discussed in Sec. IIII1 The difference 
is in that for the first three cases one has two conditions for finding the JEF solution parameters A and (3 while in the 
remaining four cases one has to satisfy one more condition. This additional constraint couples the model parameters 
Ak to the lattice spacing h. As a result, in the last four cases, for fixed Ak, one has TI solutions only at particular h, 
while in the first three cases, even for fixed Ak, one has TI solutions for any h. 
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FIG. 4: Spectrum of the lattice with a kink in the case (iii) with only A2 and A4 nonzero. Straight horizontal line at uj = \fl 
shows the lower bound of the phonon band while dots show the kink's internal modes calculated for the kink at various positions 
xo with respect to the lattice. At any position xq the kink possesses the zero-frequency Goldstone translational mode. The 
kink profiles at xq = and xo = 0.5 are shown in Fig. [3] (b) and (c), respectively. Model parameters: h = 0.8, A2 = A4 = 0.5, 
A = 1. 

Let us demonstrate this qualitative difference between two groups of models by comparison of the properties of the 
static kinks. The first group of models will be represented by the case (iii) with only Ai and A4 nonzero while from 
the four models of the second group we will choose the case (iv) with only A3 and A§ nonzero, and the case (vi) with 
only A3, A 4 , and A 5 nonzero. For the kink solutions discussed below we will always set A = 1. 

Kink in the case (iii) with only Ai and A4 nonzero. Parameters of the kink solution Eq. (|19|) with 5 = 1 for this 
case are given by Eq. (I24[) . For given model parameters h and A4 one can find the inverse kink width (3 solving 
the second equation in Eq. (|24|). The model parameter Ai must satisfy the continuity constraint given by the last 
expression in Eq. (|24f . Particular feature of this discrete <f>4 model is that it admits the TI solutions at constant Ak , 
see Fig. [3] Contrary to that, as it will be seen in the following examples, in the models (iv) to (vii), model parameters 
Ak are /i-dependent. 

Vibrational spectrum of the lattice containing a kink at different positions with respect to the lattice xq is shown in 
Fig. [4] for model parameters h — 0.8, A2 = A4 = 0.5, A = 1. The corresponding kink profiles at xo — and £0 = 0.5 
are shown in Fig. [3] (b) and (c), respectively. At any position xo the kink possesses the zero-frequency Goldstone 
translational mode. Straight horizontal line at uj = \[2 shows the lower bound of the phonon band, see Eq. (|197|) . 

Kink in the case (iv) with only A3, and A§ nonzero. Model parameters and kink parameters in this case are given 
by Eq. (|36|) . For chosen A3 (or A5) one can find A5 (or A3) from the continuity constraint [Eq. (|36[) ] and then find 
the inverse kink width (3 solving the second equation in Eq. (|36[) . Finally, one of the Eq. (|3(j[) relates the model 
parameters Ak to the lattice spacing h. In Fig. O we show (a) the model parameters Ak and kink inverse width (3 
as functions of h. In the region of lattice parameter around h — 0.5 it is possible to have two different static kink 
solutions at the same h which is illustrated in (b) and (c). In both cases h = 0.5, but model parameters Ak and the 
inverse kink width /3 are different (shown in each panel). Both kinks are stable and have a zero- frequency translational 
Goldstone mode at any position with respect to the lattice x$. 

Kink in the case (vii) with only A3, A4, and A§ nonzero. Kink and model parameters are related by Eq. (|58p . In 
Fig. [5]we plot (a) model parameters Ak as the functions of h at fixed inverse kink width (3 = 2; (b) the on-site kink at 
h = 1.3; and (c) the inter-site kink at h = 1.3. Other model parameters for (b) and (c) are A3 — 0.8297, A4 = 1.0092, 
A 5 = -0.8389, and A = 1. 

Note that in Fig. [6] (a) all Ak vary with h but, interestingly, the inverse kink width (3 is constant (= 2) in a wide 
range of lattice spacing h. In the classical discrete 4> 4 model and in the models with /i-independent parameters Ak the 
kink width usually decreases with increase in h. On the other hand, it is possible to get from Eq. (f5"8"]) one constant 
model parameter with two other model parameters Ak and kink parameter (3 being functions of h. 

In Fig. [7] one can see the spectrum of the lattice with a kink. Straight horizontal line at u = v2 shows the 
lower bound of the phonon band [see Eq. (|197[) ] while dots show the kink's internal modes calculated for the kink at 
various positions xq with respect to the lattice. At any position xq the kink possesses the zero-frequency Goldstone 
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FIG. 5: Static kink in the case (iv) with only A3 and As nonzero, (a) Model parameters Ak and kink inverse width (3 as 
functions of h. In the region of lattice parameter around h = 0.5 it is possible to have two different static kink solutions at the 
same h which is illustrated in (b) and (c). In both cases h = 0.5, but model parameters Ak and the inverse kink width /3 are 
different (shown in each panel). Both kinks are stable and have zero- frequency translational Goldstone mode at any position 
with respect to the lattice xo- 



translational mode. The kink profiles at xq = and xq = 0.5 are shown in Fig. [S] (b) and (c), respectively. Model 
parameters are: h = 1.3, A 3 = 0.8297, At = 1.0092, A 5 = -0.8389, and A = 1. 



VIII. DISCUSSION AND CONCLUSIONS 

In this paper we have introduced a rather general discrete </> 4 model of which all known models in the literature are 
special cases. We could find seven special cases when the model as given by Eq. (fTJ) supports the exact static JEF 
and hence hyperbolic kink and pulse solutions. Two of those seven cases have been analyzed in [l2[ and while 
for the remaining five cases, JEF solutions were given in Sec. Mil 

The exact solutions constructed for the considered discrete <j) A model are important for the theory of the TI lattices. 
Indeed, the JEF static solutions with an arbitrary shift along the lattice xq are the TI solutions with the zero- 
frequency Goldstone mode, i.e., solutions that are free of the Peierls-Nabarro potential. The discrete (j) 4 model of Eq. 
(fTj) corresponding to the following four particular cases 

1. case (iv) only A3 and A§ nonzero; 

2. case (v) only A2, A3, and A$ nonzero; 

3. case (vi) only A3, A4, and A§ nonzero; 

4. case (vii) only A2, A3, A4, and A$ nonzero; 

are the new TI models. Each of these models (like the other three) supports a two-dimensional set of TI static 
solutions that can be parameterized by the points of the plane (m, xq). For fixed model parameters Ak these models 
support TI solutions only for a particular lattice spacing h. Note, however, that the other three models, i.e. cases 
(i) to (iii), for fixed model parameters A2 and A4, support the TI solutions for any arbitrary value of h. However, 
when Ak are considered to be functions of h, the TI solutions in (iv) to (vii) can also be constructed for continuously 
varying h (see numerical examples of Sec. I VII B[) . In this context, it is worth noting that the TI Model 5 given by 
Eq. ([5]) and the TI Model 9 also have /i-dependent nonlinearities. 
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FIG. 6: Static kink in the case (vii) with only ^3, A4, and As nonzero, (a) Model parameters At as functions of h at fixed 
inverse kink width j3 = 2. (b) On-site kink at h — 1.3. (c) Inter-site kink at h = 1.3. Other model parameters for (b) and (c) 
are A 3 = 0.8297, Aa = 1.0092, A 5 = -0.8389, and A = 1. 



In this paper, we also showed that the general model, Eq. (|TJ) supports periodic sine and staggered sine solutions. 
Remarkably, almost all the known models (even those not supporting the JEF solutions), were found to support these 
solutions. Besides, a large number of exact, short-periodic and aperiodic static solutions admitted by Eq. ([1]) were 
obtained in Sec. IIVI While we do not have a rigorous proof, but the few examples discussed in Sec. IV CI suggest that 
very likely, the short-periodic, aperiodic as well as trigonometric solutions, in fact follow from low order algebraic 
equations. In this context, it is worth pointing out that the sine solution does not follow from the map for Model 
2 as well as the map for case (iii) when only A2 and A4 are nonzero. The factorization can also easily explain the 
appearance of the aperiodic solutions that can be regarded as the solutions obtained from different multipliers and 
linked together, as exemplified by the discussion below Eq. (|183[) . 

It is worth pointing out that none of the factorized problems discussed in Sec. IV CI contain an integration constant 
and thus they generate only particular solutions. Some of them are TI solutions, for example, the three-periodic 
solution to the Speight and Ward Model 7 derivable from Eq. (|187[) and Eq. (|188[) . while others are not, for example, 
arbitrary sequence of ±1, derivable from Eq. (|183p . 

As it was shown in Sec. IV CI in some cases the two-point map can be factorized and then the short-periodic solutions 
can also be derived from a two-point problem. We also discussed several examples in which the three-point problem 
can be reduced to a set of two lower-order finite-difference equations, and one of those equations is a two-point one 
while another is a three-point one. Based on these examples, we conjecture that all short-periodic, aperiodic as well 
as sine and staggered sine solutions should follow from lower order difference equations. 

Note that the short-periodic solutions and, more generally, the solutions derived from factorized problems very 
often do not survive the continuum limit because factorized equations usually have a different continuum limit than 
the original, non-factorized one. In this context it is worth noting that sine is not a solution of the continuum <f> field 
equation. One exception to this rule is the kink solution to the Speight and Ward Model 7 for which the reduced 
two-point problem Eq. in the continuum limit obtains a form which is equivalent to the first integral of the static 
4 field. 

Coming back to the exact JEF solutions, we emphasize that they are important because by using them one can 
construct the corresponding two-point maps from which the corresponding solutions can be obtained iteratively. 
Moreover, in some cases, the map obtained for a particular JEF solution can be transformed to the form of a general 
map from which majority of static solutions including other JEF solutions admitted by the model can be constructed. 
We conjecture that, except for those obtainable from special cases when the three-point problem can be substituted 
with a set of two equations, all other solutions can be obtained recursively from this general map. Following this way 
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FIG. 7: Spectrum of the lattice with a kink in the case (vii) with only A3, A4, and A5 nonzero. Straight horizontal line at 
u — v2 shows the lower bound of the phonon band while dots show the kink's internal modes calculated for the kink at various 
positions xo with respect to the lattice. At any position xo the kink possesses the zero-frequency Goldstone translational mode. 
The kink profiles at xo — and xo = 0.5 are shown in Fig. [6](b) and (c), respectively. Model parameters: h = 1.3, A3 = 0.8297, 
A 4 = 1.0092, A 5 = -0.8389, and A = 1. 



we could construct the map Eq. (|158|) from which any static solution of Eq. (fTJ) with only A2 and A4 nonzero can be 
constructed (again, except for the solutions that result from specially factorized three-point problems). On the other 
hand, for cases (iv) to (vii), while one can obtain a map from a JEF solution, so far we are unable to obtain a general 
map. 

In Sec. IVII Al we provided numerical evidence that the Speight and Ward Model 7 does not support TI static 
solutions other than those derivable from reduced lower-order algebraic problems as discussed in Sec. IV CI In 
addition to the well-known TI kink solution we have found the TI sin and staggered sin solutions to this model. We 
believe that in the general (non-factorized) formulation, the static Speight and Ward model is not integrable and a 
two-point map that includes the integration constant as a parameter cannot be constructed for this model. 

Based on the results of the present study one can separate the TI models into two classes. In the first class belong 
the models that support a two-dimensional space of TI static solutions. These solutions, if they are derivable from a 
two-point nonlinear map, can be parameterized by the points of the plane (C, <fio), where C is the integration constant 
that can vary continuously within certain range and 4>q is the initial value of the map that can also vary continuously. 
Alternatively, if the JEF solutions are known, then the solutions can be parameterized by the points of the plane 
(m,Xo) so that m plays the role of the integration constant C while variation of xq plays the role of </)rj, and results 
in the shift of the solution along the lattice. The second class is formed by the models that admit TI static solutions 
with an arbitrary shift along the lattice (controlled by either xq or cf>o) but corresponding solutions do not include the 
integration constant as a parameter. 

The TI models in the first class have been investigated in [l2|, UH, UH and four more TI cubic nonlinearities belonging 
to this class are found in the present work. In particular, TI I models are cases (i) to (vii) and Model 2. In this 
context it is worth noting that while a universal two-point map is known for Model 2 (for arbitrary 7 and S), no JEF 
or any other analytical solutions are known so far which can be characterized by C and </V On the other hand, no 
universal two-point map is known for cases (iv) to (vii). 

It is likely that the Speight and Ward model 7 [5| belongs to the second class of TI models because it supports 
the well-known TI kink and the TI sine solutions derived in the present work [see Eq. (|52")l and Eq. ([57]) ] but 
these solutions are derived from reduced equations, as shown in Sec. IV CI The reduced equations do not contain 
the integration constant. On the other hand, e.g., the five-periodic solution derived from the non-factorized model 
possesses the Peierls-Nabarro potential, as shown in Sec. IVII Al 

At this stage, it may be worthwhile classifying the various known models. As noted above, cases (i) to (vii) and 
Models 2 and 9 (and hence Models 3 to 6) are TI I models. On the other hand, Models 7, 8 and 10 arc in general TI 
II models even though in the special cases ofev = /3 = Oora; = 0,/3 = —1/2, Model 8 becomes a TI I model. 

Note that while a TI I model can also be a TI II model, the converse is obviously not true. For example, except for 
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case (ii) when only A4 is nonzero, all other TI I models also admit sine and staggered sine solutions. Besides, Model 2, 
Model 8 (in case 7 = 0) and Model 10 (in case a>i = 0:3) also satisfy the four-periodic TI solution (..., a, b, —a, —b, ...). 
At the special value of A = 2, such a TI solution is admitted by Model 9 (at j3 = 0) as well as by Models where only 
A3, A4, A§ or only A2, A3, A4, A§ are nonzero. 

It is thus clear that the general Model as given by Eq. ([1]) is only a TI II model even though in few special cases, it 
could be elevated to a TI I model. There is one model, however, which is not a TI model with respect to any known 
solution. We have in mind the model where only A\ is nonzero, i.e. the <f> model with standard discretization. For 
example, while it admits the short-period solutions (..., a, 0, —a, ...), (..., a, 0, —a, 0, ...), (..., a, a, —a, —a, ...), neither of 
them possesses the Goldstone mode. 

Before closing, we spell out some of the open problems. 

1. For the cases (iv) to (vii), can one obtain a unified general two-point map from which all solutions, including the 
JEF solutions can be derived? Note that, at the moment by starting from sn, cn, dn solutions, one can obtain 
three different maps from which only the respective solution can be obtained. 

2. For Model 2, while a general two-point map is known, to date no analytic solution is known which is characterized 
by the two parameters C and <f>$. Can one find few such analytic solutions? 

3. Can one rigorously show that all short period, sine and staggered sine solutions for any TI I or TI II model, 
follow from the lower order equations? 

4. Can one rigorously prove that Model 7 of Speight and Ward is only a TI II and not a TI I model? 

5. There is a belief that all TI models (at least TI I models) must have some conserved quantity. Unfortunately, 
for cases (iii) to (vii) (which are all TI I models) no such conserved quantity is known at present. Can one find 
such a quantity or disprove the conjecture? 

6. While it has been demonstrated that no discrete model can simultaneously have conservation of P\ and energy 
E, it is not known whether one can have a model where both P2 and E can be simultaneously conserved. The 
obvious guess would be no. It would be nice to prove or disprove this conjecture. 

7. Can one find particular TI solutions for the discrete models that do not belong to the TI I class, i.e., finding 
the isolated TI solutions to the discrete models that are not considered by many researchers as the TI models. 
In the present study we have given several examples of such solutions, for instance, four-periodic TI solution 
(..., a, b, —a, —b, ...) and the TI trigonometric solutions given in Sec. IIII CI It can be so that isolated TI solutions 
exist for many discrete models. Many of the isolated TI solutions result from factorized static problems and 
thus, finding various factorizations of the original static problem can be a method for their derivation. 

Finally, the results obtained in this paper are easily extended to the case of the general nonlinear Schrodinger 
equation [2l|. We hope to address these issues in a forthcoming publication. 
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IX. APPENDIX 

In this Appendix, using the various short-period solutions given in Sec. IV, we spell out the solutions which 
are admitted by some of the Models discussed in the paper. We shall only mention the allowed solutions with the 
lowest period. Needless to say that the corresponding solutions with arbitrarily large period or aperiodic solutions 
(if they exist) will also be valid in that case. For example, in case only A2 is nonzero, as shown below, the solution 
<f> n = (..., a, a, —a, ...) is allowed provided A = 2 and o 2 = 1. It is then clear from the discussion in Sec. IV that in that 
case a solution of the form <j) n — (..., a, —a, (a, —a p times), —a, ...) with period 2p + 1 (p > 1) as well as an aperiodic 
solution with any number of "a" and "-a" kept at random but with the constraint that at most two "a" or two "-a" 
are always together, is also an exact solution. 
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A. Model 1: Only A\ nonzero, Aih 2 = A 

This is the case of the classical discretization of <fi A term. In this case the admitted solutions are 

(i) (fi n — (..., a, —a, ...) in case a 2 — ^j-^- Thus such a solution is valid if A > 4 or if A < 0. 

(ii) (f> n = (..., a, —a, 0, ...) in case a 2 = ^j^-- Thus such a solution is valid if A > 3 or if A < 0. 

(hi) 4> n — (..., a, a, —a, —a, ...) and (fi n = (..., a, 0, —a, 0, ...) in case a 2 = Thus such a solution is valid if A > 2 

or if A < 0. 

(iv) 4> n = (..., a, a, 0, —a, —a, 0, ...) in case a 2 = ^j^- Thus such a solution is valid if A > 1 or if A < 0. 

B. Model 3: Only A 2 nonzero, A 2 h 2 = A 
In this case the admitted solutions are 

(i) 4> n = (..., a, —a, ...) in case a 2 = ^j^-- Thus such a solution is valid if < A < 4. 

(ii) (f> n = (..., a, —a, 0, ...) in case a 2 = ^j^-- Thus such a solution is valid if < A < 3. 

(hi) <fi n = (..., a, a, —a, —a, ...) and (fi n — (..., a, 0, —a, 0, ...) in case A = 2 while a is cm?/ real number. 

(iv) (f> n = (..., a, 6, —a, —6, ...) in case A = 2 while a, b both are any arbitrary real numbers. 

(v) <fi n — (..., a, a, 0, —a, —a, 0, ...) in case a 2 = ^-i. Thus such a solution is valid if A > 1 or if A < 0. 

(vi) In addition, there are a large number of solutions which are admitted in case A = 2 and a 2 = 1. These allowed 
solutions are 

(fi n = (..., a, a, — a, ...) , (..., a, a, a, — a, ...) , (..., a, 0, — a, a, ...) , (..., a, a, a, — a, — a, ...) , 

(..., a, — a, a, — a, 0, ...) , (..., a, a, — a, — a, 0, ...) , (..., a, a, 0, — a, 0, ...) , (..., a, 0, — a, a, a, ...) , 

(..., a, a, a, 0, —a, —a, ...) , (..., a, a, 0, — a, 0, a, ...) , (..., a, —a, a, a, —a, 0, ...) , (..., a, a, 0, — a, — a, 0, a, ...) , 

(..., a, 0, —a, a, 0, —a, 0, ...) , (..., a, — a, 0, a, a, 0—, a, ...) , (..., a, a, —a, a, a, —a, 0, ...) , (..., a, 0, — a, a, a, a, —a, ...) , 

(..., a, — a, 0, a, a, a, 0, —a, ...) , (..., a, a, 0, —a, —a, 0, a, —a, 0, ...) , 

(..., a, a, 0, — a, — a, 0, a, — a, 0, a, ...) , (..., a, —a, 0, a, a, a, 0, — a, a, a, — a, ...) . (205) 

C. Model 6: Only A4 nonzero, A 4 h 2 = A 
In this case the admitted solutions are 

(i) (fi n — (..., a, —a, ...) in case a 2 = ^j^-- Thus such a solution is valid if either A > 4 or A < 0. 

(ii) (fi n = (..., a, a, —a, —a, ...) in case a 2 = ^j^-- Thus such a solution is valid if < A < 2. 

(hi) (fin = (••■, a, 1/a, a, 1/a, ...) and c/>„ = (..., a, 0, —a, 0, ...) in case A = 2 while a is any real number, 
(iv) (fin = (■••, a, —a, 0, ...) in case A = 3 while a is any real number. 

In addition, there are a large number of solutions which are admitted in case A = 1 and a 2 = 1. These allowed 
solutions are 

(fin = (..., a, a, a, -a, -a, ...) , (..., a, a, -a, -a,0, ...) , (..., a, a, 0, -a, -a,0, ...) , (..., a, a, a, 0, -a, -a, ...) , 

(..., a, a, 0, — a, — a, 0, a, ...) . (206) 

D. Case (iii): Only A 2 ,A 4 nonzero, (A 2 + A 4 )h 2 = A 

In this case the admitted solutions are 

(i) (fin — (••■, a, —a, ...) in case h 2 a 2 (A4 — A 2 ) = A — 4. 

(ii) (fi n = (..., a, 1/a, a, 1/a, ...) in case h 2 A 4 = 2 , h 2 A 2 = A — 2 while a is any real number. 

(iii) (fin = (•••, a, —a, 0, ...) in case a 2 h 2 A 2 = 2(3 — A). 

(iv) (fin — (■••, a, a, —a, —a, ...) in case a 2 h 2 A4 = 2 — A. 

(v) (fin — (■■■, a, 0, —a, 0, ...) in case A = 2 while A 2 , A4 as well as a are arbitrary real numbers. 

(vi) (fin = (..., a, a, —a, —a, 0, ...) in case a 2 h 2 A 2 = 2A — 2. 

(vii) c/>„ = (...,a,a,0, -a,0, ...) , (..., a, 0, -a, a, 0, -a, 0, ...) , 
(..., a, a, 0, —a, —a, 0, a, —a, 0, ...) , in case A = 2, h 2 a 2 A 2 = 2. 

(viii) (fin = (..., a, a, a, —a, —a, ...) , (..., a, a, —a, —a, 0, ...) , (..., a, a, 0, —a, —a, 0, a, ...) , in case a 2 = 1, h 2 A 2 = 2A — 
2, h 2 a 2 At = 2 - A. 
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(ix) 4> n — (..., a, a, —a, ...) , (..., a, —a, a, —a, 0, ...) , (..., a, —a, a, a, —a, 0, ...) , in case 

2 8-3A - _2A(3-A) u2 A A(2 - A) 



Thus this solution is valid provided < A < 8/3. 



E. Case (iv): Only A 3 ,A 5 nonzero, (A 3 + A 5 )h 2 = A 



In this case the admitted solutions are 

(i) 4> n — (..., a, —a, ...) in case /i 2 a 2 (A 3 — A 5 ) = A — 4. 

(ii) <j> n = (...,a,0, ...) in case a 2 h 2 A 5 = 2, A = 2. 

(iii) </>„ = (..., a, —a, 0, ...) in case a 2 h 2 A 3 = 2(A — 3). 

(iv) 0„ = (..., a, a, 0, ...) in case a 2 — 1, /i 2 A 3 = A — 1, /i 2 ^ = 1. ( v) <f> n = (..., a, 0, —a, 0, ...) in case A = 2 and a is 
any real number. 

(vi) (f> n = (..., a, a, —a, —a, ...) in case a 2 h 2 A 3 = A — 2. 

(vii) 0„ = (..., a, a, 0, —a, 0, ...) in case A = 2, a 2 /i 2 A 3 = 2. 

(viii) 4> n — (..., a, a, 0, —a, —a, 0, ...) in case a 2 h 2 A 3 = 2A — 2. 

(ix) 4> n = (..., a, a, a, 0, -a, -a,0, ...) in case a 2 = 1, /i 2 yl 3 = 2(A - 1), h 2 A 5 = 2 - A. 

(x) (j) n — (..., a, 0, —a, a, 0, —a, 0, ...) in case A = 2, a 2 h 2 A 3 = —2. 

(xi) <f> n = (..., a, —a, 0, —a, a, 0, ...) in case 



o 2(A-2) , 2/ , A(A-3) A 



a ' 3 = iv-ir ' 5 = x^2 • (208) 



Thus this solution is valid provided either A > 2 or A < 0. 

(xii) <j) n = (..., a, —a, a, —a, 0, ...) in case 

2 _3A_8 2 _ 2A(A-3) 2 _ A(A - 2) 
a ft A 3 - 3A _ g , hA 5 -^—^. (209) 

Thus this solution is valid provided either A>8/3orA<0. 

(xiii) <t> n = (..., a, a, —a, ...) , (..., a, a, a, —a, ...) , (..., a, a, a, —a, —a, ...) , in case a 2 — 1, /i 2 A 3 = A — 2, h 2 a 2 A 5 = 2. 

(xiv) In addition, there are a large number of solutions which are admitted in case A 3 = A 5 , A = 4, a 2 = 1. These 
allowed solutions are 

<f> n = (...,a, -a,0,a, ...), (..., a, a, -a, 0, -a, ...) , (..., a, -a, a, a, -a, 0, ...) , (..., a, a, -a, 0, -a, a, ...) , 

(..., a, 0, — a, a, a, a, — a, ...) , (..., a, a, — a, a, a, — a, 0, a, ...) . (210) 



F. SW Model 7: Ai = f , A 2 = A 3 = f , A 4 = A 5 = 0, A 6 = A 
In this case the admitted solutions are 

(i) </>„ = (..., a, —a, ...) in case a 2 = £L^zfl . Thus such a solution is valid if either A > 4 or A < 0. 

(ii) 4> n = (..., a, —a, 0, ...) in case a 2 = 6 ^ A ~ 3 ^ . Thus such a solution is valid if either A > 3 or A < 0. 

(iii) 4> n = (..., a, a, —a, —a, ...) in case a 2 = 9 ^^" 2 - ) . Thus such a solution is valid if either A > 2 or A < 0. 

(iv) tp n — (..., a, 0, —a, 0, ...) in case a 2 = 9( - 2 ^ 2 - ) . Thus such a solution is valid if cither A > 2 or A < 0. 

(v) cf> n = (..., a, a, 0, —a, —a, 0, ...) in case a 2 = 18 ^ A . Thus such a solution is valid if either A > 1 or A < 0. 

(vi) (j) n = (..., a, a, 0, —a, 0, ...), (..., a, a, 0, —a, 0, a, ...) in case a 2 = 1, A = 18/7. 

(vii) 4> n = (..., a, a, 0, ...), (..., a, a, —a, —a, 0, ...) in case a 2 = 3/2, A = 12. 

(viii) (j> n — (..., a, a, —a, ...), (..., a, a, a, —a, —a, ...) in case a 2 = 1, A = 9/2. 

(ix) <f> n = (..., a, a, 0, —a, —a, 0, a, —a, 0, ...) in case a 2 = 6/5, A = 15/4. 

(x) In addition, there arc a large number of solutions which are admitted in case A = 6 and a 2 = 3. These allowed 
solutions are 



4> n = (.-, a, 0, ...) , (...,a,-a,a,-a,0,...), (..., a, 0, a, 0, -a, ...) , (..., a, -a, 0, -a, a, 0, ...) , 
(..., a, 0, a, 0, a, —a, ...),(..., a, 0, —a, a, 0, —a, 0, ...) . 



(211) 
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G. Model 2 with 8 = 0, 7 = 1/4: Only A 5 , A 6 nonzero with A 5 h 2 = A 6 /i 2 = A/2 

This is same as the Model 4 due to Kevrekidis. In this case the admitted solutions are 

(i) (j> n = (..., a, —a, ...) in case a 2 = ^j^-- Thus such a solution is valid if < A < 4. 

(ii) (j> n = (..., a, —a, 0, ...) in case a 2 = Thus such a solution is valid if < A < 3. 

(hi) (j) n = (..., a, a, —a, —a, ...) and <p n — (..., a, 0, —a, 0, ...) in case A = 2 while a is any real number. 

(iv) 4> n = (.-., a, o, —a, —6, ...) in case A = 2 while a, b both are any arbitrary real numbers. 

(v) 4> n = (•••) a i a 7 0, —a, —a, 0, ...) in case a 2 = 4 ^ A A ^ . Thus such a solution is valid if A > 1 or if A < 0. 

(vi) In addition, there are a large number of solutions which are admitted in case A = 2 and a 2 = 1. These allowed 
solutions are 

(j) n = (..., a, a, -a,...), (..., a, a, a, -a, ...) , (..., a, 0, -a, a, ...) , (..., a, a, a, -a, -a, ...) , 

(..., a, — a, a, — a, 0, ...) , (..., a, a, — a, — a, 0, ...) , (..., a, a, 0, — a, 0, ...) , (..., a, 0, — a, a, a, ...) , 

(..., a, a, a, 0, —a, —a, ...) , (..., a, a, 0, — a, 0, a, ...) , (..., a, —a, a, a, —a, 0, ...) , (..., a, a, 0, — a, — a, 0, a, ...) , 

(..., a, 0, —a, a, 0, —a, 0, ...) , (..., a, — a, 0, a, a, 0—, a, ...) , (..., a, a, —a, a, a, —a, 0, ...) , (..., a, 0, — a, a, a, a, —a, ...) , 

(..., a, — a, 0, a, a, a, 0, —a, ...) , (..., a, a, 0, —a, —a, 0, a, —a, 0, ...) , 

(..., a, a, 0, —a, —a, 0, a, —a, 0, a, ...),(..., a, —a, 0, a, a, a, 0, —a, a, a, —a, ...) , (212) 



H. Model 2 with 8 = 1/4, 7 = 0: Only A U A 3 , A 4 nonzero with A x h 2 = A 4 h 2 = A :i h 2 /2 = A/4 
In this case the admitted solutions are 

(i) 4> n — (..., a, —a, ...) in case a 2 = ^j^- Thus such a solution is valid if cither A > 4 or A < 0. 

(ii) <f>„ = (..., a, —a, 0, ...) in case a 2 = 2 ^ A r 3 ^ . Thus such a solution is valid if either A > 3 or A < 0. 

(iii) <j) n = (..., a, b, —a, —6, ...) in case (a 2 + b 2 ) = 2 ^ 2 ^ = 2. 

(iv) 0„ = (..., a, a, a, —a, ...) in case a 2 = 1, A = 4. 

(v) 0„ = (..., a, a, 0, —a, 0, ...) in case a 2 = 4/3, A = 3. 

(vi) <f> n = (..., a, a, 0, —a, —a, 0, a, ...) in case a 2 = 1, A = 2. 

(vii) In addition, there are a large number of solutions which are admitted in case a 2 = 2 ^ A A ~ 1 ^ . These allowed 
solutions are 

4> n = (..., a, a, -a,...), (..., a, a, -a, -a, ...) , (..., a, 0, -a, 0, ...) , (..., a, a, 0, -a, -a, 0, ...) . (213) 



I. Model 2 with arbitrary 5,7: Ai = A 4 = A 3 /2 = 5X,A 5 = ^ 6 = 27A, A 2 = (1 — 4<S — 47) A 

This is the full Model 2. Note that we have already obtained solutions in three special cases. In this general case, 
many more solutions are admitted. The admitted solutions are 

(i) 4> n = (..., a, —a, ...) in case a 2 = ^Efjx ■ Thus such a solution is valid if either A > 4 or A < 0. 

(ii) 4> n = (..., a, 0, ...) in case 

a 2 = ^-^, 2 7 Aa 2 = l. (214) 
oA 



(iii) 4> n = (..., a, —a, 0, ...) in case 



(iv) <f> n — (..., a, a, —a, ...) in case 



(v) <p n — (■••, a, a, 0, ...) in case 



2 3A-4 A-2 
° =-A-' <5= 2(3A-4)- (216) 



fl2 =27A = I^- ^ 



(vi) (j) n — (..., a, a, a, 0, ...) in case 

(vii) <p n — (..., a, a, a, —a, ...) in case 

(viii) <f>„ = (..., a, b, —a, —6, ...) in case 

(ix) <f> n = (..., a, a, —a, —a, ...) in case 

(x) 4> n = (..., a, 0, —a, 0, ...) in case 

(xi) 4> n = (..., a, —a, 0, a, ...) in case 



a 2 = 1 , <*> = -!-, 7 = 

2 ' ' 2A 



a 2 = 1, 5 = 0, A = 2. 



/ 2 ,2n A - 2 

(a 2 + b 2 ) 



SA 



2 A - 2 



2 A-2 
a 2 = 



<5A 



a 2 = 7r , A = 3 , 85 + 2 7 = 1 . 
07 



(xii) (f> n — (..., a, a, —a, —a, ...) in case 

A-2 



a 2 = 1 , (5 = 



2A 

(xiii) </>„ = (..., a, —a, 0, a, ...) in case 



a 2 = -!- , A = 3 , 8(5 + 2 7 = 1 . 

6 7 



(xiv) </>„ = (..., 0, a, 0, a, a, ...) in case 



■ A+2 r A-2 1 



a =— r~> = o/a 1 1 7 = 



A ' 2(A + 2)' ' 2(A + 2)' 

(xv) <ft n — (..., a, —a, a, —a, 0, ...) in case 

fl2 = l7T' (8<5-l)(A-2) = 2 7 (A-4). 

(xvi) 4> n = (..., a, 0, a, 0, —a, ...) in case 

2 3(2A-3) A-2 1 



A ' 3(2A-3)' ' 6(2A 
(xvii) <fi n — (..., a, a, —a, —a, 0, ...) in case 



,2 



(A-2) 1 



2<5A 1 - 4 7 

(xviii) </>„ = (..., a, a, 0, —a, 0, ...) in case 



(xix) <f> n — (..., a, —a, 0, —a, 0, ...) in case 



2 (2A-1) . A- 2 
a =" * *=7^7a TT ' 7 



A ' 2(2A-1)' ' 2(2A-1) 

(xx) (fi n = (..., a, a, 0, a, —a, ...) in case 

» 3 4 A = 3 - <->-w 

(xxi) </>„ = (..., a, a, 0, —a, —a, 0, ...) in case 

2 2(A-1) 



a 



(1 + 26 - 4 7 )A ' 

(xxii) <fi n = (..., a, —a, 0, —a, a, 0, ...) in case 

a 2 =2^A> (85-l) = 2 7 (2A-7). 

(xxiii) (fi n = (..., a, a, a, 0, —a, —a, ...) in case 

a 2 = l, 7 = 0, *=^- 

(xxiv) (fi n — (..., a, a, 0, —a, 0, a, ...) in case 

2 A-2 A-2 

a = 1> 7 = < * , o = — - — . 

4A A 

(xxv) (fin — (■•■, a, a, 0, a, —a, 0, ...) in case 

2 _ (6A + 7) 4A - 7 5 

° ~ 5A ' 2(6A + 7) ' 7 ~ 2(6A + 7) 

(xxvi) </>„ = (..., a, a, —a, 0, —a, a, ...) in case 

a 2 = l, A = l, 7 =i, <5=^- 



(xxvii) 0„ = (..., a, 0, a, 0, a, —a, ...) in case 



a 2 = 3, A = 3, 7=^j, ^=^- 

(xxviii) (fin = (•••, a, 0, a, 0, a, a, ...) in case 

a 2 = 1 , A = 2 , 7=1, 5 = 0. 

(xxix) (fin — (..., a, —a, a, a, 0, a, ...) and (..., a, a, 0, a, a, 0, —a, ...) in case 

(xxx) (fin — (..., a, a, 0, —a, —a, 0, a, ...) in case 

,2 - - - A - 2 



a 2 = 1 , 5-27 



2A 

(xxxi) 0„ = (..., a, 0, —a, a, 0, —a, 0, ...) in case 

« 2 = ^, 25(5-3A) = (27-l)(A-2). 
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(xxxii) <ft n — (..., a, a, 0, —a, —a, a, —a, ...) in case 



.. (7A - 12) A - 2 _ 7(A - 2) 



a ~ A ' S ~ 2(7 A - 12) ' 7 ~ 4(7A - 12) ' (244) 

(xxxiii) In case 

a 2 = 1, (5 = 7 = 0, A = 2, (245) 

then the following solutions are allowed 

4> n = (..., a, 0, —a, a, a, ...),(..., a, 0, —a, a, a, a, —a, ...), (..., a, — a, 0, a, a, a, 0, —a, ...), (..., a, 0, a, 0, a, a, a, — a, ...) , 
(..., a, —a, 0, a, 0, a, a, 0, —a, ...),(..., a, —a, 0, a, 0, a, a, a, 0, —a, ...),(..., a, —a, 0, a, a, a, 0, —a, a, o, —a, ...) . (246) 

(xxxiv) tfin — (..., a, a, 0, a, 0, —a, —a, ...) in case 

a 2 = 2, <5 = 0, 7 =-L. (247) 

(xxxv) 4>n — (■••, a, a, —a, 0, —a, 0, —a, ...) in case 

9 5-2A „ 1 , , 

° =— <* = °' ^2(5^21)- « 



(xxxvi) 4> n — (..., a, —a, a, a, a, 0, —a, ...) in case 



a 2 = 1, (5 = 0, 7 = A = 2. (249) 



(xxxvii) 4> n — (..., a, a, 0, —a, —a, 0, a, —a, 0, ...) in case 



« 2 = Tn 57 — 5~^r i 6(l-37) = (3A-l)(l-3<y-37). (250) 



A(l - 35 - 3 7 ) 
(xxxviii) </>„ = (..., a, a, 0, —a, —a, 0, a, —a, 0, a, ) in case 



.2 1 x 7(A - 2) _ 1 



(xxxiv) In case 



1 ' 5 = 18A"^ ' 7 =M- (251) 



, 3A-4 A- 2 , N 

° = ^' 5 = 7 =2(3Ar-l)' ( 252 ) 



then the following solutions are allowed 

4> n = (..., a, 0, —a, a, ...),(..., a, —a, a, a, —a, 0, ...),(..., a, —a, 0, a, a, 0, —a, ...),(..., a, a, —a, a, a, —a, 0, ...) . (253) 



J. Case (v): Only A 3 , A4, A 5 nonzero, (A s + A 4 + A 5 )h 2 = A 

In this case the admitted solutions are 

(i) 4> n — (..., o, —a, ...) in case h 2 a 2 (A 3 + A 4 — A 5 ) = A - 4. 

(ii) 4> n = (..., a, 0, ...) in case a 2 h 2 A 5 = 2 , A = 2. 
(hi) 4> n — (..., a, —a, 0, ...) in case a 2 h 2 A 3 = 2(A — 3). 

(iv) 4> n = (..., a, a, -a, ...) in case a 2 h 2 (A 3 - A A ) = (A - 2),a 2 h 2 (H + F - J) = A - 4. 

(v) 0„ = (..., a, a, 0, ...) in case a 2 /i 2 A 3 = 2(A - l),a 2 h 2 A 5 = 2. 

(vi) = (..., a, 0, —a, 0, ...) in case A = 2 and a is any real number. 

(vii) <f) n = (..., a, a, —a, —a, ...) in case a 2 h 2 A 3 = A — 2. 

(viii) 0„ = (..., a, a, a, 0, ...) in case a 2 = 1, h 2 A 3 = 2(A - 1) , h 2 A 5 = 2, h 2 A 4 = -A. 

(ix) cj) n = (..., a, -a, 0, a, ...) in case a 2 = 1, h 2 A 3 = 2(A - 3), h 2 A A = 4 - A, h 2 A 5 = 2. 

(x) <j) n = (..., a, 6, —a, —6, ...) in case A3 = A4, A = 2 and a, b are arbitrary real numbers. 
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xi ) 4 


>n = (••• 


a, a, 0, —a 


,0, ...) 


in case A 


= z, h a A3 = 2. 






■h — ( 


- (J, (Z 7 a, 


-a n 

a, u, 


1 in r*fiQP 

. . t 111 LclOC 


h 2 a 2 ( An + /t, — AA - 


= A — i 


xiii) 


4>n = (• 


., a, a, a, — 


a, —a, 


...) in case 


a 2 = l,/i 2 (A - A) 


= A- 


xiv) 


4>n = (• 


., 0, a, 0, a 


a, ...) 


in case A 


= A, A = 2,fr 2 a 2 A 


= 2. 


xv) 


K = (.. 


, a, 0, a, 0, 


-a,... 


in case A 


= 2, A 3 = -A 5 ,h 2 A 4 


= 2, ft 


xvi) 


= (• 


., a, a, —a, 


-a,0, 


...) in case 


h 2 a 2 A 4 = A, h 2 a 2 A 3 


= 2(A 



xvii) 0„ = (..., a, —a, 0, —a, a, ...) in case 

3A-8 , 9 . 2A(A- 



A 



h 2 A 3 = 



3) 



3A-8 



h 2 A, = 



A(A-4) 
3A-8 



2. 



^ 2 A = 



2A 



3A-8 



xviii) </>„ = (..., a, a, 0, —a, —a, 0, ...) in case a 2 h 2 A 3 = 2(A — 1). 



xix) (j> n = (..., a, — a, 0, — a, a, 0, 

xx) (/>„ = (•.., a, a, a, 0, —a, —a, ...) in case a 

2 



) in case a 2 h 2 A 3 = 2(A - 3), h 2 a 2 A 5 = 2. 



1, A 3 = -A 5 , h 2 A 3 = 2(A - 1), /i 2 A = A. 

xxi) </>„ = (..., a, a, 0, —a, 0, a, ...) in case a 2 = 1, A = 2, /i 2 A = 2. 

xxii) 0„ = (..., a, —a, a, a, —a, 0, ...) in case 



a 2 = 



5A- 16 



A3 = A. 



h 2 A 3 = 



2A(A-3) 
5A- 16 



^ 2 A - A(A ~ 4) 
" A4 _ 5A - 16 



(xxiii) cp n — (..., a, a, 0, —a, —a, 0, a, ...) in case a 2 = 1, h 2 A 3 = 2(A — 1). 

(xxiv) (f> n — (..., a, 0, —a, a, 0, —a, 0, ...) in case A = 2, a 2 h 2 A 3 = —2. 

(xxv) (p n = (..., a, a, 0, —a, —a, a, —a, ...) in case A = 2, h 2 a 2 A 3 = 2 , h 2 a 2 (A 4 — A) = — 4. 

(xxvi) 0„ = (..., a, a, 0, a, 0, -a, —a, ...) in case A = 2, a 2 = 3, h 2 A 3 = h 2 A = /i 2 A = 2/3. 

(xxvii) 0„ = (..., a, —a, 0, a, 0, a, a, a, 0, —a, ...) in case A = 2, h 2 a 2 A 3 = 2, A = 2A + A- 



(254) 



(255) 
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